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Abstract 

We examine the potential of using colliders to distinguish models with parity {Z2) stabilized dark matter 
(DM) from models in which the DM is stabilized by other symmetries, taking the latter to be a Z3 symmetry 
for illustration. The key observation is that a heavier mother particle charged under a Z3 stabilization 
symmetry can decay into one or two DM particles along with Standard Model (SM) particles. This can be 
contrasted with the decay of a mother particle charged under a parity symmetry; typically, only one DM 



particle appears in the decay chain. The larXiv: 1003.08991 studied the distributions of visible invariant mass 
from the decay of a single such mother particle in order to highlight the resulting distinctive signatures of 
Z3 symmetry versus parity symmetry stabilized dark matter candidates. We now describe a complementary 
study which focuses on decay chains of the two mother particles which are necessarily present in these events. 
We also include in our analysis the missing energy/momentum in the event. For the symmetry stabilized 
mothers, the resulting inclusive final state can have two, three or four DM particles. In contrast, models 
with Z2 symmetry can have only two. We show that the shapes and edges of the distribution of MT2-type 
variables, along with ratio of the visible momentum/energy on the two sides of the event, are powerful 
in distinguishing these different scenarios. Finally we conclude by outlining future work which focuses 
on reducing combinatoric ambiguities from reconstructing multi-jet events. Increasing the reconstruction 
efficiency can allow better reconstruction of events with two or three dark matter candidates in the final 
state. 



1 Introduction 



A stable weakly interacting massive particle (WIMP) - with a mass also of order the weak scale 
- is a well-motivated candidate for dark matter (DM) in the Universe since it approximately has 
the correct relic density upon thermal freeze-out [1]. Such a particle also often arises in extensions 
of the Standard Model (SM), especially those motivated by solutions to the Planck- weak hierarchy 
problem of the SM. Finally, if the WIMP is a part of an extension of the SM, then it is likely to 
have (weak) interactions with SM particles. Hence, the WIMP paradigm can be tested via non- 
gravitational methods, for example, direct /indirect detection of cosmic DM or production of the 
DM at colliders. The latter is our interest here. 

The collider searches of the DM paradigm typically involve producing a heavier particle charged 
under the same symmetry which stabilizes the DM. Such a "mother" particle must decay to SM 
particles and DM, manifesting as missing energy, along with a SM final state. Reconstructing the 
decay chains leading to such events will enable us to determine the masses of the DM and the mother 
particles; it is advantageous to do this in a model- independent manner. A tremendous amount of 
effort has been put-in into such a research program, especially at the CERN Large Hadron Collider 
(LHC), see, for example [SlElllElEllZlElEllIOllIIllIllliaillllSlliaiia 
See also [23] for a review. Most of this work still assumes that a Z2/parity symmetry stabilizes 
the DM (henceforth called Z2 models). This is partly because there are several models with such 
a symmetry \25\ [26} [271 l28l 129] , However, while Z2 might be the simplest possibility for such a 
symmetry, it is by no means the only one [301 1331 ESI IM] ■ So, in previous work [311 [3S], the program 
of instead determining the DM stabilization symmetry from collider data was started. We focused 
on how to distinguish a DM stabilization symmetry (henceforth called a Z^ model) from a Z2 
modej^. We have focused on Z3 models for simplicity and definiteness. However, we emphasize our 
techniques in this and in previous papers [311 [36] can be generalized to distinguish most other DM 
stabilization symmetries from parity symmetry stabilized DM. 

The basic idea behind distinguishing Z3 from Z2 models is that a 

• single mother charged under a Z3 symmetry is allowed (based simply on the symmetry) to 
decay into one or two DM candidates. 

This is to be contrasted with the fact that mother particles charged under a Z2 symmetry have 
only one DM candidate in the final stat^. As discussed in reference [36], decays of a single Z3- 
charged mother particle generate a "double edge" in the invariant mass distribution of the visible 

^See reference [3S] for how calculation of relic density and indirect detection might be modified in a Z3 model vs. 
a Z2 model. Also see reference [31] for a collider study of non-Z2 models, but with decays of mother particles into 
DM occurring outside the detector vs. the case of such decays taking place inside the detector studied here and in 
reference [36j . 

In both Z3 and Z2 models, a mother can decay into three DM particles, but we will not consider this possibility 
here for simplicity and since it is expected to be phase-space suppressed compared to the other possibilities. 
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(SM) particles. This is with the condition that the intermediate particles in the decay chains are 
off-shell and that the decay chains with one and two DM contain identical SM particles. For the 
case of on-shell intermediate particles, it was shown in [36] that this invariant mass distribution 
has a "cusp" for certain decay topology (with two DM) of a Z3 mother particle. In all, ^6] focused 
on new features in observables from a single decay chain only. 

In this work we consider the total inclusive event in order to glean even more information, 
recalling that there must be two such mother particles present. For example, consider the case 
where there is only one visible (SM) particle in the decay chain of a mother particle. Constructing 
the invariant mass of the visible particle of this decay chain, as per the analysis of |36j . is not 
very useful for the purpose of reconstructing the mass of the mother particle: one might have to 
resort to including information about the invisible particle(s) in the same decay chain. Since we 
can only measure the total missing transverse momentum in the event which is shared between 
invisible particles from two mothers, we must use measurements from both sides. An option is to 
use "MT2"-type observables/variables [2l[3lll[5l[6l[71[8l[9l[ini[IIl[I21[I3[I2[l5l[l6] Another case 
where one of the analyses of reference [3S] (based on single mother decay) might not work is when 
the visible/SM particles in the decay chains with one and two DM (of course for Z3 model) are not 
identical (even if they are more than one). Thus, one does not obtain a double edge for the case of 
intermediate particles in the decay chains being off-shell. 

With the above motivations in mind, in this paper, 

• we develop techniques for distinguishing models from Z2 models using information from 
both mother decays and the missing (in addition to visible) energy/momentum in an event. 

We especially study the above cases where the techniques of reference [36] might not work - in this 
sense, our present work is complementary to that of reference [36] . 

• We show that shapes and edges of these Mt2 distributions, along with the ratio of visible 
momentum/energy on the two sides of the event, act as powerful discriminants between Z3 
and Z2 models (including the case of a neutrino, i.e., massZess invisible particle - in addition 
to DM, in the final state for Z2 models). 

An outline of this paper is as follows: In the next section, we begin with a review of Mt2 variable 
in Z2 models. We present some important formulas such as the location of maximum value of the 
Mt2 distribution, and discuss some interesting features such as a "kink" in the maximum Mt2 as 
a function of the "trial" DM mass in the case where more than one visible particle is involved in 
each decay chain. In Sec. [3] we move on to the Mt2 variable in Z3 models. We define three different 
event types arising in Z3 models based on the number of dark matter particles in the final state 
and provide their corresponding theoretical predictions of the maximum Mt2- In particular, we 
discuss the conditions to have a kink in maximum A'It2 as a function of trial DM mass for the cases 
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where there exist more than two DM particles in the final state: such situations do not arise in Z2 
models. We further show some simulation results for M'T2 distributions for the new types of events 
in Z3 models and discuss some notable features to be used for distinguishing Z3 models from Zi 
models. 

In the following two sections, we provide detailed applications of our results of the previous two 
sections for distinguishing Z3 from Z3 models. First we consider the easier case where the decay 
chain with one DM contains visible/SM particle(s) which are not identical to the ones in the two 
DM decay chain mentioned above (in Z3 models). Based on the theoretical considerations given in 
Sees. [2]and[3l we provide ways of distinguishing Z3 models from Z2 models, as well as measuring 
the mother and DM masses. We do it for both the case of one visible/SM particle in the two decay 
chains and more than one visible/SM particle case. In the next section we deal with the case where 
one DM and two DM decay chains contain identical visible/SM particle(s), and discuss additional 
techniques required in this case to distinguish Z2 and Z3 models. In all of the examples above, we 
make the simplifying assumption that the intermediate particles in the decay chain are ojff -shell, 
i.e., they are heavier than the mass of their mother particle. In sectional we briefly mention some 
aspects of the case of intermediate particles being on-shell. We next conclude and show details of 
some of the calculations in the Appendix. 

2 A Review of Mt2 for Z2 Models 

For simplicity, in this paper 

• we consider pair-production of a single type of mother particle which is charged under the 
DM stabilization symmetry. 

We also assume that the total transverse momentum of the two mother particles produced in 
an event is zero, for example, we neglect any initial/final state radiation. In Z2 models, each 
such mother decays into SM/visible particle(s) and one DM/invisible particljf . Furthermore, it is 
assumed 

• we know which visible particle(s) originate from which decay chain0 For example, if the 
pair-produced mother particles are boosted sufficiently, their decay products are likely to be 
collimated so that the visible particles coming from the same decay chain are detected in the 
same hemisphere in the collider. For alternate methods of determining the correct assignment 
of visible particles to the two decay chains, see reference [37] and section 16. 3[ 



We assume that there is only one type of DM particle in this (and similarly the Z3) model so that the invisible/DM 
particles in each decay chain are identical. 

^Of course, for Z3 models, we do not know which decay chain emits one or two DM particles. 
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The Mt2 variable [21 [3] is a generalization of the transverse masjfl to this case. Specifically, for 
each event, it is defined to be a minimization of the maximum of the two transverse masses in each 
decay chain under the constraint that the sum of all the transverse momenta of the visible and 
invisible particles vanishes [21 [3] : 

Mr2= min [max(M^^\ M^^Hl (1) 

1.(1), 1,(2) -(1) -(2) „ L L ^ J J 

where p^^*^ denote the vector sum of visible transverse momenta and p^^ denote the transverse 
momentum of the invisible particle in the ith decay chain (i = 1, 2): the minimization is performed 
over the latter momenta. Mj. is the usual transverse mass: 

[M^^f = [mf^f + m2 + 2 (sf^ - p^^'^ • p«) (2) 

where E'^^*^ and m^^*^ are (respectively) the transverse energy and transverse mass formed by 
all visible particles belonging to the same decay chain. The variables with a tilde represent the 
corresponding quantities formed by the invisible particle in the same decay chain. Note that the 
mass of the invisible particle rh should be regarded as a unknown/free parameter because we are 
not aware of it in advance, and henceforth we call it "trial" DM mass. In this sense Mt2 should be 
considered to be a function of the trial DM mass m, and its maximum value among many events 
(which will be used extensively in the following) is defined as 

M^^''{m) = max [MT2(m)] ■ (3) 
many events 

Obviously M^2^{m) is also a function of the trial DM mass (see App.|A|for details). An important 
result to be noted is that if there are a sufficient number of events and the actual DM mass is 
substituted into m, then the above-given M^2^ becomes the actual mass of the pair-produced 
mother particles [21 E] : 

M^^"" {m^mDM)^M (4) 

where M and indicate the true masses of mother and DM, respectively. 

In order to see how this Mt2 analysis is applied to realistic situations, we first take the case 
where there exists a single visible/SM particle in each decay chain, and then move on to the case 
where there exists more than one visible/SM particle in each decay chain. A similar analysis can 
be done for the mixed case, i.e., one visible particle on one side and more than one on the other. 

2.1 One Visible/SM Particle in Each Decay Chain 

In this case the upper edge in Mt2 distribution is obtained by "balanced" [21 El S] solution (see 
App. [A] for details). 



Mt2 -M^2 -V -^2 +V 4M2 ^"^ 



^Of course, the usual transverse mass assumes only a single mother particle. 
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Figure 1: Theoretical expectation of M^2^ versus the trial mass m for Z2 events. The masses of 
mother and DM particles are 400 GeV and 100 GeV, respectively. The left panel shows the case 
where there exists only a single visible particle per chain. The right panel shows the case where 
there exists more than one visible particle per decay chain. In both panels, the solid black curve 
represents the overall/net upper edge, M^*^^. In the right panel, the dotted straight line which 
extends into the right-hand part of the solid line is the M^^^ for the unbalanced solution, whereas 
the dashed line which extends into the left-hand part of the solid line is that for the balanced 
solution. 

Here (and henceforth) 

• we assume that all visible particles are massless for simplicity. 

As mentioned earlier, the above-given upper edge is a function of the trial DM mass m and one 
can see that it reduces to the true mother mass M with rh equal to the true DM mass moM- The 
left panel of Fig. [T] shows the above theoretical prediction for the location of maximum Mt2 for 
Z2 models as a function of the trial mass m. We used 400 GeV and 100 GeV as mother and DM 
particle masses. As expected from Eq. ([5]), the curve in the figure "smoothly" increases with fh 
(cf. the following section), and that the value corresponding to m = moM (here 100 GeV) is 

the same as the true mother particle mass (here 400 GeV: see the black dotted lines). 

2.2 More than One Visible/SM Particle in Each Decay Chain 

Once there exists more than one visible particle per decay chain, another type of solution to Mt2, 
denoted by "unbalanced" [2l[3lll], arises. If 

• we assume that the intermediate particles in the decay chains are Oj^-shell, i.e., heavier than 
their mother particles 
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(as we will for most of this paper), then the balanced solution in this case is still given by Eq. ([5]), 
and the unbalanced solution is as followa^: 



rriDM + m. 



(6) 



Hence the overall upper edge in the Mj'2 distribution is determined by a "competition" between 
balanced and unbalanced solutions: 



max Mrp^ , Mr 



'■T2 



4M2 



for 771 > m DM 



for rh < rriDM- 



(7) 



Note that the M^^^ shows different functional behaviors depending on the relative size of the trial 
DM mass to the true DM mass. As a result, M^^^ is no longer smoothly increasing with m in 
contrast to the case with one visible/SM per decay chain. Instead, there arises a "kink" at the 
location of the actual DM mass, with the corresponding M^2^ being the actual mother mass 
This is illustrated in the right panel of Fig. [T] where the upper edges for the two possible types of 
solutions ^^^^^'^^^"■^ QT^d M^^^'^'^^ are shown by a straight line (i.e., the dotted line, which extends 
into the right-hand part of the solid line) and a dashed curve (which extends into the left-hand part 
of the solid curve), respectively. The upper edge in the Mt2 distribution is given by the larger of 
these two values, i.e., the black solid curve, and shows a kink (indicated by the black dotted lines), 
at m = 100 GeV and M^2^^ = 400 GeV (as expected). 



3 Mt2 for ^3 Models 

To begin with, we would like to reiterate some of the relevant features of Z3 symmetry in order to 
avoid any possible confusion later. Under the Z3 symmetry, a particle/field transforms as 



)exp 



(8) 



where g = (i.e., neutral) or +2. Suppose the lightest of the Z3 charged particles (labeled <^o) 
has charge g' = +1 (similar argument can be made for charge q - +2). Clearly, its anti-particle (<^o) 
has (a different) charge q--l (which is equivalent to q- +2) and has the same mass as (/>o. Then, 
solely based on symmetry considerations, all other (heavier) Z3-charged particles can decay into 
this lightest Z3-charged particle or its anti-particle (in addition to Z3-neutral particles, including 
SM particles). To be explicit, a heavier Z3-charged particle with charge q - +1 can decay into 
either (single) (f)Q or two 0o's (and Z3-neutral particles). Taking the CP conjugate of the preceding 



Of course, in general, one can find the expressions for both the balanced and the unbalanced solutions in the case 
of oft-shell intermediate particles [3]. 

'A similar kink also appears for the case of one visible particle in each decay chain if the total transverse momentum 
of the two mother particles is non-zero, for example, in the presence of initial/final state radiation (see [5]), but (as 
mentioned earlier) we neglect this possibility for simplicity. 
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statement, we see that a heavier Za-charged particle with the other type of charge, namely q - 2, 
is allowed to decay into two 0o's or single (pQ. Of course, (po or cannot decay and thus is the 
(single) DM candidate in this theory. We will denote this DM particle and its anti-particle by DM 



and DM, respectively, although we do not make this distinction in the text since DM and anti-DM 



particles are still degenerated 

According to the above- given argument, for Z3 models, each mother particle can emit either 
one or two DM particles so that there exist two, three, or four DM particles in the final state (for 
pair-production of mothers) while there are only two DM particles for Z2 models. We therefore 
expect richer structures in the Mt2 distribution for Z3 models. Here 

• we take as an ansatz only a single DM particle in each decay chain for the sole purpose of 
defining My, even if there could be two DM particles in either or both of the two decay chains. 

We do so for the following two reasons. First, in the real collider experiment, there is (a priori) 
no clear information on the number of invisible particles involved in the decay process of interest so 
that each individual decay chain with only one DM is a natural (starting) assumption. Moreover, 
one can naturally expect (and we will show) that decay events from Z3 models will show different 
features in the Mt2 analysis compared with those for Z2 models. Therefore, starting with a Z2 
assumption and deriving a "contradiction" in the Mt2 analysis, we can distinguish Z^ models 
from Z2 ones (which is our primary goal here). In this context, we call such an analysis imposing 
one-DM-per-chain assumption the "naive" Mt2 analysis. 

For a more systematic consideration let us define the three different events having different 
numbers of (a single type of) DM particles as E2, -E3, and £^4-type, respectively, i.e., each subscript 
on E simply implies the total number of DM particles in the final state: see Fig. [21 where SMi^ 2 

le actual number of particles in the state) in the 
Here the red dashed lines denote any particles 
charged under dark matter stabilization symmetry (in this section Z^ symmetry) while the black 
solid lines/arrows denote any visible/SM particles. One should note that E4 type events represent 
the case with 2 DM particles in each decay chain. Also, both decay chains (with one and two DM) 
might not exist for a specific mother so that all three types of events might not occur. Like in Z2 
models, we start with the case with one visible/SM particle in each decay chain, and we consider 
the case with more than one visible/SM particle in each decay chain in the following subsection. 
While doing so, we see how the Mt2 analysis applied to Z3 models contrasts with Z2 models. 



denote the visible/SM final states (irrespective of t 
decay chains with one and two DM (respectively]^ 



*Of course, which of the two particles is denoted anti-DM is a matter of convention. Also, as a corollary, the DM 
particle should be Dirac fermion or complex scalar in a Z3 model. 

^These two SM final states might not or might be identical: we will return to these two possibilities in the next 
two sections (respectively). 
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i) type ii) type iii) E^ type 



Figure 2: The three types of events m Z3 models, based on the total number of DM in the event. 
"M" denotes the mother particle. Each SM final state can have more than one particle, and the 
subscripts 1 and 2 on SM denote the visible/SM final states in the decay chains with one and two 
DM, respectively. Note that, based simply on the Z^, symmetry, if a mother decays into DM, then 
the same mother decays into two aniz-DM in the other decay chain. Since DM and anti-DM have 
same mass and are not detected, we neglect this distinction between the two henceforth (we already 
did so thus far). For simplicity, we will also henceforth not differentiate between SM and SM or 
between M and M. 



3.1 One Visible/SM Particle in Each Decay Chain 

In this case the upper edge in the (naive) Mt2 distribution is determined only by the balanced 
solution like Z2 models, and the analytic expressions for the three type events are given as follows 
(see App. [X]for details): 



Mt2,E2 - Mt2,E, -y -^2 +V UP ^ ^ ' 



As a reminder, the events with E2^ and E4 represent events with two, three and four dark matter 
candidates. Note that Eq. ([9]) has the same form as Eq. ([5]) in Z2 models because E2 type events 
also contain two DM particles in the final state (just like Z2 models) whereas the other two types 
of events do not appear in Z2 models so that the corresponding Eqs. (jlOp and (jlip (and similar 
ones later) are new/do not appear in previous literature. Substituting rh - munj in Eq. ([9]) gives 
the true mother mass M for the value of M^2^ (as expected), but the other two equations give a 
combination of the true mother and DM masses rather than the true mother mass. Actually, this is 
not surprising because we have used the naive Mt2 variable for the S3 4-type events in Z3 models, 
whereas the actual physics is different from the physics under which our Mt2 variable is defined. 
For example, for £"3 type events there is an asymmetry between the final states of the two decay 
chains, which is caused by adding one more DM to either of the two decay chains. For £'4 type 
events, even though the two decay chains have symmetric final states, the "effective" DM mass is 
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Figure 3: Theoretical expectation of M^2^ versus the trial mass rh for Z3 model. The masses 
of mother and DM particles are 400 GeV and 100 GeV, respectively. The left panel shows the 
case where there exists only a single visible particle per chain. The black, red, and blue curves 
are showing the corresponding MJ^2^ values to £"2, -Ba, and type events over m, respectively. 
The right panel shows the case where there exists more than one visible particle per decay chain. 
The overall upper edges, MJ^2^ for E2 and i?4 type events, are given by the solid black and blue 
curves, whereas the balanced and unbalanced solutions are denoted by the dashed and dotted curves 
(respectively) which merge into the solid curves on the right (left)-hand part. The corresponding 
plot for £"3 type events can be found in the next figure. 



twice the true DM mass so that the true mother mass (for the value of M^^^) is in fact obtained 
by setting m = 2m om instead as clearly seen from Eq. (Ilip . 

All of the theoretical predictions mentioned above are demonstrated in the left panel of Fig. O 
Again, we used 400 GeV and 100 GeV as mother and DM masses. The black solid, the red dashed, 
and the blue dot-dashed curves represent the theoretical expectations of M^^2^ for Ej,, and £4 
type events, respectively (the curve for £'2-type events is of course the same as in right-hand side 
of Fig. [T]). As discussed above, the M^2^ for E2 type events corresponds to the true mass of the 
mother particle (here 400 GeV) with the trial DM mass equal to the true DM mass (here 100 GeV) 
whereas -B4 type events do it for twice the DM mass (here 200 GeV), as shown by the dotted/black 
lines. 

In addition, there are a couple of features to be noted; there is no kink arising in the M^*" 
curves for £3^4 -type events just like the case of a single visible particle per decay chain in Z2 models 
(or an £'2-type event in Z3 models). Second for any given m, the values form a hierarchy of 

M^^E,>M^^%^>M^^^^. (12) 
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3.2 More than One Visible/SM Particle in Each Decay Chain 

Once more visible particle(s) are added in each decay chain, one could naturahy expect that a kink 
appears like in Z2 models. The reason is that, just like for Z2 models, the maximum unbalanced 
solutions take part in determining the overall upper edge in the Mt2 distribution together with the 
balanced solutions. It turns out, however, that this expectation is true only for E2 and E/^ type 
events which we discuss to begin with. Again, assuming the intermediate particles are ojff -shell the 
maximum values of the balanced solutions for E2 and £'4 type events are simply given (as for the 
one visible particle case) by Eqs. ([9]) and (jlip . respectively, and those of the unbalanced solutions 
are given as follows: 

M^g"^^""^"' = M-m£,M + fh for £'2 (13) 

^max,«n6a« ^ M-2mDM + m ioT E^. (14) 

Here Eq. (|13|) is of exactly the same form as Eq. ^ due to the similarity between the decay 
structures for Z2 models and E2 type events while Eq. (jl4p for £'4-type events is relevant only for 
Z3 models, i.e., it is not present in Z2 models. More quantitatively, the above-given two equations 
differ by munj for any given m because one more DM particle is emitted in both decay chains for 
Ei type events compared with E2 type events (see Eqs. and (pl|) in App. [A|) . 

As mentioned for Z2 models, the maximum Mt2 values are given by the larger of the balanced 
and unbalanced solutions: for E2 type events as in Eq. ([7]), and for E4 type events by 

AI - 2m DM + ITT' for m > 2m dm 

(15) 



Mmax [ n rmax,bal n ^max,unbal 

T2,E4 = maxlM^2,E4 ' ^T2,E4 



(M2-4m2„.,)2 /(M2-4m2„.,)2 _ 
^ 4A/" + V iM-^ + ^ fo'^ ^ ^ ^rUDM- 



Note that there is also a kink for i?4-type events as seen from Eq. (|15|) . but aX ffi- 2mDM (i-e., not 
tudm), because the effective DM mass for £'4 is given by 2mDM- 

The right panel of Fig. [3] illustrates the above theoretical considerations for E2 and £4 type 
events (of course the curve for E2-tYpe events is the same as the right-hand side of Fig. [1]). As 
before, the two straight dotted lines which extend into the right-hand parts of the solid lines indicate 
the maximum values of the unbalanced solutions while the two dashed curves which extend into 
left-hand parts of the solid curve indicate the maximum values of the balanced solutions. The 
actual upper edge in the Mt2 distribution for any m is given by the black (for E2 type events) or 
blue (for E4 type events) solid curves. Identifying the location of the kink in £2-type events and 
its corresponding M^2^ enables us to determine the masses of mother and DM particles separately 
(just like in Z2 models). The figure also shows the kink for ii'4-type events, but located at m = 200 
GeV (i.e., 2mDM, as expected) and M^2^ - 400 GeV. This observation can be used as a cross-check 
for the determination of M and moM based on £'2-type events (again, this feature is new in Z3 
models relative to Z2). 
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Figure 4: Theoretical expectation of M^2^ versus the trial mass rh for type events. The mass of 
mother particle is 400 GeV for both cases, but the masses of DM particle to be used are 150 GeV 
and 100 GeV for the left panel and the right panel, respectively. For both cases, the black solid 
lines give the maximum of Mt2, whereas the dashed curves give the maximum for the balanced 
solution. The maximum for the unbalanced solution coincides with the solid line on the left panel, 
whereas on the right panel, it is given by the dotted straight line (which extends into the right-hand 
part of the solid line). 



On the other hand, as far as type events (which are absent in Z2 models) are concerned, 
whether or not there exists a kink depends on the mass hierarchy between mother and DM particles 
(see App.|B]for details). Again, assuming Oj^-shell intermediate particles the maximum balanced 
solution is simply given by Eq. (fTOl) (just like the case of one visible particle per decay chain), 
whereas the maximum unbalanced solution has the same form as that for E2 type events because 
one of the two decay chains still emits a single DM particle in the final state (see Eqs. ()64p and 
(PI) in App.m. 



max, unfeai 
T2,E3 



M - rriuM + rh 



for E^ 



(16) 



If the ratio of the DM mass to the mother mass is larger than {\/3 - l)/2, it turns out that 
the maximum unbalanced solution given in Eq. ()16p is always bigger than the maximum balanced 
solution given in Eq. (fTUj) so that 

V3- 



r2,ij3 =maxlMy2,E3 ^ ^T2,Ez 



M - m DM + "rh 



for > 



1 



and for all m. (17) 



M 2 

The left panel of Fig. [5] clearly confirms our expectation (based on above equation) that there 
occurs no kink in the upper edge of Mt2 as a function of the trial DM mass, i.e., the upper edge 
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in the Mt2 distribution is always determined by the unbalanced solution (black solid line), not by 
the balanced solution (red dashed curve). Here we adopted M - 400 GeV and m^M - 150 GeV, 
and thus the ratio between them is obviously larger than (Vs - l)/2. 

On the other hand, once the ratio of DM to mother mass is smaller than (^/3 - l)/2, the 
competition between the balanced and the unbalanced solutions results in 



I ^-max 



max, 6a/ , rraa.x.unbal 

M - m OM + "fn 



and tor — — < 



4AP ^ '"' 



for rh>m' 
^ for m<m' 



M 



(18) 



where 



m 



(M - niDM) (^(M2-m2,^,)(M2-4m2,j,,) - M{M - msA/)) 



2M{M-mDM)-y/{AP' 



m 



DM 



(19) 



We see that there is a kink at rh-m'. Here m' is not simply the true DM mass but it is given by 
a combination of the true mother and DM masses (in fact, it is smaller than the true DM mass), 
which is clearly different from that in Z2 models. The functional behavior of M^2^ for this case is 
shown in the right panel of Fig. 31 Here we took rrioM - 100 GeV which makes the ratio smaller 
than (Vs - l)/2. As before, the maximum Mt2 for the balanced and unbalanced solutions are 
shown by the dashed and dotted curves (which extend into the black solid curve to the RHS and 
LHS). The final maximum Mt2 is given by the larger of these two solutions (black solid curve) 
which clearly shows a kink at a value of m which is different from the actual DM mass ttldm - 100 
GeV (shown by the vertical black dotted line) as expected based on above discussion. Of course, we 
can still evaluate the masses of mother and DM particles (using £^3-type events only) by obtaining 
M^2^ and m' from the above Mt2 analysis, substituting them into Eqs. (jlSp and ()19p . and solving 
those two equations about M and moM- 

Next, let us investigate the hierarchy among the three M^2^ values for E2, E^, and i?4 type 
events. Although a bit more complicated than the one visible particle case, it is nonetheless 
straightforward to derive this hierarchy based on above equations, We have a following hierarchy 
among the M^^^ values for the three types (cf. the one visible particle case shown in Eq. (jl2p ): 



for rh > niDM 
for 771 < rriDM- 



(20) 



13 



E2 Events 



£3 Events 
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Mt2 (GeV) 

Figure 5: The Mt2 distributions for E2 (top left) and E^, (top right) type events and E2 + 1' events 
(bottom). The mother and the DM particle masses are 400 GeV and 150 GeV, respectively, and 
the trial DM mass (?fi) used is 6 GeV. 

3.3 Shapes of Mt2 Distributions 

Before closing the present section, let us examine the shape of the Mt2 distributions for Z3 cases. 
For this purpose we simulated events using MadGraph/MadEvent [3H]. Here and in sections 15.1.11 
and l5.1.2t we make the following assumptions (mostly for simplicity): (a) effects of spin correlations 
are neglected in both production and decay of mother particles, i.e., we assumed scalar particles 
only in our simulations,; (b) the beam is proton-proton with 14 TeV total energy in the center 
of mass frame (motivated by the LHC parameters); (c) the non-colored (scalar) mother particles 
are pair-produced via s-channel exchange (of another scalar particle); and finally (d) only the 
relative values of the number of events (on vertical scale) are meaningful, i.e. only the shape of the 
distributions is the robust aspect of our analysis. 

For all the simulations in this section, the masses of mother and DM particles used in the (toy) 
model are 400 GeV and 150 GeV, respectively, we took a single visible particle per decay chain for 
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simplicity and chose a small trial DM mass rh (here 6 GeV) for the purpose of illustration. The 
upper left panel of Fig. [5] demonstrates the Mt2 distribution for E2 type events. Obviously, this is 
similar to the Mt2 distribution for events in Z2 models [21 [3l 3] , which is not surprising because 
E2 type events also have only two dark matter particles in the final state like Z2. The (naive) 
Mt2 distribution for E3, type events illustrated in the upper right panel of Fig. [5] has two notable 
features. As expected from the analytic expressions given before, first of all, the location of the 
upper edge is clearly lower than that for E2 type events (or for Z2 models). Secondly, the shape of 
Mt2 distribution for E^ type events shows a long tail near the upper edge compared with E2 type 
events (which have relatively sharp upper edge): this is because more physical constraints (e.g., 
rapidity) between decay products should be satisfied in £^3 type events in order that they form 
a kinematic configuration to give maximum Mt2, thereby reducing the corresponding number of 
events near ^1^2"^. 

This feature of a (relatively) long tail is also true for a special case in Z2 models with three 
invisible particles in the full event one of which is a (massless) SM neutrino, i.e., with one DM in 
one decay chain, but one DM and a neutrino in the other (henceforth we call it E2 + 1' events!^. The 
bottom panel of Fig. [5] demonstrates the Mt2 distribution for E2 + 1' events, assuming intermediate 
particles are off-shell lj. At the upper edge it shows a tail which is longer than that for E2 type 



events (or neutrino-less events in Z2 models), but which is not as long as that for £^3 type event j^^l. 
Such a shape is not surprising because the additional invisible particle (neutrino) m. E2 + y vs. 
E2-type events is massless so that E2 + y events can be understood as a transitional type between 
E2 and E^ type events. However, it is crucial to note that the location of the upper edge for the 
E2 + V events is the same as that for usual (i.e., neutrino-less) events in Z2 models with only a 
single massive invisible particle per decay chain, because the effective DM mass in former decay 
chains is also moM- We can therefore distinguish the E2 + 1' events from the -E3-type events in Z3 
models by observing the location of the upper edge: the latter events will have a smaller edge. 

4 Applications: Non-Identical Visible Particles in the Two Decay 
Chains 

Next, we apply the theoretical observations on the Mt2 technique for Z2 and Z^ models, which are 
described in the previous sections, for distinguishing models from Z2 ones in some specific cases. 
Like in the previous sections, we assume (for simplicity) pair-produced (same) mother particles and 
that all visible particles are massless and use the naive Mt2 analysis for all events. As mentioned 
earlier, in Z3 models each mother particle can decay into either one or two dark matter particles 



Mt2 distributions for such events have been studied in [3] and in |24) . 

^^As for the events with three (or more) DM, we assume here that there is a single (massive) invisible particle in 



each decay chain for the purpose of defining Mr 

^^This fea 
decay chain 



^^This feature of a long tail in E2 + i^-type events is valid even for the case of more than one visible particle in each 
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along with visible/SM particle(s). Here, 



• we consider a mother particle in a Z3 model for which both these decay chains (with one and 
two DM, respectively) exist. 

Similarly, in Z2 models the mother particles can only decay into a single dark matter particle along 
with visible/SM particle(s). Here, we assume two such decay chains for a Z2-mother which have the 
same visible final states as the above two Zs-mother decay chains (respectively). The idea behind 
this choice is that such a Z2 model could easily fake a Z3 one (at least based on the identity of 
the visible states). This motivates us to distinguish these two types of DM stabilization symmetry 
using the Mt2 variable. 

For later convenience, we divide the discussion into two cases based on whether or not the 
visible state in the decay chain with one DM is identical to the one in the decay chain with two 
DM in the Z3 model. We begin with the case where the SM final states in the two decay chains 
are not identical (this includes the case of partial overlap between these final states). Following 
the notation of previous section, let us denote SMi^ 2 to be these SM final state particles - whether 
they consist of one or more SM particles - in the two decay chains. And, for the Z3 model, assume 
that SMi comes with one DM particle and SM2 is associated with two DM particles!^ We thus 
have three distinct (based simply on identity of visible states) types of events in both the Z3 and Z2 
models, denoted by SMu, i.e., SMi on each side/from each mother and similarly SM12, 22- Clearly, 
for Z3 models, these three types of events correspond (respectively) to i?2,3,4-type events mentioned 
in the previous section and shown in Fig. [2] (with SMi ^ SM2). Hence, we can apply the formulas 
for the theoretical predictions of M^2^ for E2, E3, and type events derived in the previous 
section to the SMu^ 12, 22 events. 

4.1 One Visible/SM Particle in Each Decay Chain 

As a further subcase, we assume that SMi^2 consist of only one particle. Clearly, the upper edges 
of the Mt2 variable for SMn, SM12, and SM22 are given by Eqs. ©, (HHD, and ([II]), respectively. 
As is obviously seen from the left panel of Fig. [3] or equivalently Eq. p2p , the location of the upper 
edge for SM12 or SM22 is lower (for all trial DM mass m: cf. the case of more than one visible 
particles below) than that for SMn (for the same mother and DM masses). In contrast, in the Z2 
model we have the same expression for M^^^ given by Eq. ([5]) for all of SMn, SM12, and SM22 
because they all involve two DM particles in the final state. Thus, 

• different edges for the SMn^ 12, 22 events (in particular, larger for SMii) can be evidence for 



Z3 models, i.e., they provide discrimination between Z2 and Z3 model^ |. 

^■^Since a Z2 model does not allow two DM particles in each decay chain, SMi and SM2 are both emitted with only 
one DM in the final state in this model. 

^■^unless multiple mother particles in the Z2 models decay into the identical final state. 
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We can be further quantitative: 

• for Z3 models we can measure the masses of mother and DM separately as follows (in spite 
of absence of kink in left panel of Fig. [3|). 

Note that the theoretical formulas for M^2^ in Eq. ([9]) through Eq. (jlip - considered as a function 
of trial mass fh - have a structure of \/C + VC + fh? where C is a constant. So, the idea is to choose 
an arbitrary trial mass, then calculate the corresponding M^2^ from the experimental data and 
thus determine the above-defined C . Also, in our specific case where visible particles are assumed 
massless, the theoretical formula for each C is written only in terms of the mother mass and the 
DM mass: 

{M'-mlM? ((M— 

4M2 4(M^2%)2 ^ ' 

4M2 4(Mf2^^^)2 ^ ' 

Solving the above-given two equations, we obtain both mother and DM masses 0: 

M = —^^aCe2-Cez) (23) 

ruDM = —^^J{ACe2 - Ce^){Ce2 - Cez) (24) 

o\/Le2 

This situation is somewhat like the double-edge signal for single mother decay studied in reference 
[36], where it was again possible to obtain mother and DM masses from the two edges in invariant 
mass distribution of visible/SM final state. 

On the other hand, for Z2 models, we obtain only a combination of mother and DM masses 
from the (single) measurement of M™2^^, given by 

V ^Ml_ ^ \y T2 — )_ = c^ 25 

and thus it is not possible to determine mother and DM masses separately. 

4.2 More than One Visible/SM Particle in Each Decay Chain 

In this case, there is more interesting behavior of M^2^ than the case of one visible particle per 
chain 0. The upper edges of Mt2 for SMu^ 12, 22 events are now obviously given by Eqs. ([7]), (jl7p . 



^^The upper edge from SM22 provides redundant information, but of course can be a cross-check. 

^^In the SM12 events, the endpoints of the visible invariant mass distributions for the two sides of the event/decay 
chains will be different in Z3 models, i.e., in _E3-type events (vs. being the same in Z2 models), already providing a 
discrimination between the two types of models [36]. However, developing another technique for distinguishing Z3 
from Z2 models based on Mt2 can still be useful. 
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(|T8j) and p5|) . respectively. As discussed in Eq. (pO|) . the upper edge for SM12 is the same as that 
for SMii (cf. one visible particle case above) for m > rriDAi, but is lower for m < itidm than that 
for SMii. And, the upper edge for SM22 events is always lower than SMu. This fact enables us to 
distinguish Z3 models from Z2 ones because in Z2 models the upper edges for SMu^ 12, 22 coincide 
for all m (just like the case with one visible particle in each decay chain). 

Moreover, there occurs a kink in the upper edge of Mt2 as a function of the trial DM mass as 
discussed in the previous sections. Because of the existence of this kink structure, SMu itself is 
sufficient for the purpose of determining mother and DM masses (again, unlike one visible case): 
the trial mass which gives rise to a kink is the true DM mass and its corresponding M^2^ is the 
true mother mass. (Of course this is how one can measure the masses of mother and DM particles 
separately even in Z2 models.) Such a direct measurement of mother and DM masses leads us to 

• a prediction (cf. one visible case) on the location of the upper edges for the other two types of 
events, namely SM12 and SM22, a confirmation of which can provide evidence for Z3 symmetry 
as underlying phvsicJ^. 

For SM12, i.e., E^-type events, actually, there are more interesting aspects of the kink structure 
in M^2^ due to the dependence on the ratio of DM and mother masses: as discussed in Sec. [3l the 
critical ratio is given by 

M 2 ^ ' 

The kink is present only when the ratio of DM and mother masses is less than the above-given 
critical ratio. In this case, the kink location can be predicted by substituting mother and DM 



masses measured from the kink in SMu into Eq. (jl9p so that it can provide a further verification 
for Z3 symmetry 

Finally, note that another way to distinguish Z3 from Z2 models in this case was discussed in 
[36]. The idea is to use SM12, i.e., £'3-type, events in Z3 models, where edges in invariant mass 
distributions of visible particles on each side are different, i.e., (M - tuom) and (M - 2 mDAi) (vs. 
the two edges being the same for Z2 models). 

4.3 Signal Fakes by an (Effective) 2nd DM Particle 

In the two previous sections we have focused on decay processes with a single type of DM particle in 
the final state (for both Z^ and Z2 models). The crucial observation for the sake of discriminating Z3 
from Z2 models is that Z3 models have more event-topologies (i.e., E2, E3, and -E4 type events with 
different upper edges in Mt2 distributions) than the case of Z2 (which has a single upper edge). 



^'^We can also predict (and then verify) location of kink in SM22 events. Alternatively, we can use kink in SM22, 
i.e., i54-type events to determine the mother and DM masses and then make predictions. 

Alternatively, this kink can be used to determine mother and DM masses, which are then used to predict 



edges/kinks in other events. 
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regardless of the number of visible/SM particles in each decay chain. In turn, this contrasting 
feature is due to the different possibilities in each decay chain in Z3 models, i.e., presence of one or 
two DM (unlike only one DM in Z2 case). 

However, Z2 models can also acquire such different possibilities for decay chains (and thus fake 

signals in the Mt2 analysis) if we assume that there is a second DM (obviously .^2-odd) particle 
(with larger mass) denoted by DM' into which the Z2-mother can decay, i.e., there are actually 
two (absolutely) stable DM particles in a Z2 model |40j . Clearly, even with only one DM in a Z2 
model, a similar effect can arise from a mother decaying into an (Z2-odd) on-shell color /electrically 
neutral particle which decays (into DM particle and SM, possibly visible), but outside the detector 
(i.e., there exists a Z2-odd particle - other than the DM ~ which is stable and invisible as far as 
the detector is concerned). Another related possibility is that there is a Z2-odd (on-shell) neutral 
particle which decays inside the detector, but invisibly, i.e., into DM and invisible SM, for example, 
neutrino. A classic example of the last type is found in supersymmetry where sneutrino decaying 
into neutrino and lightest neutralino (which is assumed to be the lightest supersymmetric particle, 
i.e., DM). Even in the latter two cases, there is "effectively" (i.e., as far as the collider analysis 
is concerned) a second "DM" and so we will denote it also by DM'. In particular, in the last 
case mentioned above, i.e., even if there is an (on-shell) neutral particle decaying invisibly inside 
detector, the theoretical prediction of the Mt2 variable is the same as with a DM' of the same mass 
as this neutral mother particle. 

Here, we note that the reference [TU] has studied such (asymmetric) events using an Mt2 type 
analysis, in particular, variants of the usual Mt2 variable have been developed. As before, we will 
instead apply the naive/usual Mt2 variable, i.e., assume (again, just for the purpose of constructing 
Mt2) that there is a single and same DM in both decay chains. 

In more detail, the above case in Z2 models gives rise to "Za-faking" signals is as follows. As 
before, consider pair production of a single mother such that decay chains with DM and DM' are 
both allowed. Consequently, we obtain three distinct decay topologies for the full event: two DM, 
one DM and one DM' and two DM'. We will denote these three types of events by £'21 -^31 ^''^^ 
E'^ since they obviously resemble (and thus can fake) £'2, £3, and £4 type events being found in 
Z^ models, respectively. In particular, we can expect three different upper edges for Mt2 in E'2, 
-E3, and £'4-type event a ^^1. Explicitly, the maximum balanced Mt2 solutions for them (for both the 
cases with one visible particle per decay chain and more than one visible particle per decay chain) 

'^''An extreme case is when DM' is massless, for example, SM neutrino. However, in this case, as mentioned in 
section 13.31 the upper edges for £3 (denoted hy E2 + v event in section 13. 3p and £'4-type events will be same as for 
£'2-type events. Thus this case can be easily distinguished from Zz models. 
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are simply given as follows: 



^''^T2,i?^ V 4^j2 +V 4M2 ^zyj 

and the maximum unbalanced Mt2 solutions (only for the case with more than one visible particle 
per decay chain) are given as follows: 

^max,nnW ^ M-niDM + rh ioT E'^ (30) 

^max,nnW ^ M - TUdM + m ioT E'^ (31) 



M^^T = M-m^A^ + m for ^4'. (32) 



Here m'jyj^{> uidm) denotes the mass of the second DM-like particle. Again, all three types 
actually contain two DM/DM- like particles, i.e., the subscripts on E' do not imply the number of 
DM particles in a full decay chain but rather indicate the respective topologies in Z3 models which 
they fake. Note that if we set m^^^ = 2mDHf, then the above edges are exactly the ones in a Z3 
model (see Eqs. ([9|) through ([TT]) and Eqs. (fT3]) . ([H]) . and (fT6|) from previous sections. This feature 
is as expected since for E-^^ 4-type events (in models) which are at the edge of the respective 
Mt2 distributions, the two DM from the same mother are collinear so that their invariant mass is 
2moMi i-e., the decay chain with two DM effectively has single DM of this mass as far as Mr2-edge 
is concerned. 

Despite the fact that such Z2 events with a second DM-like particle can introduce three decay 
topologies, we can still differentiate Z3 and Z2 models. However, the strategies to be applied depend 
on the number /identity of visible particles in each decay chain. In this section, we consider the case 
where the visible particles in the decay chain with DM (denoted by SMi, following the notation 
used earlier) are different than the visible particles (denoted by SM2) in decay chain with DM' 
in the Z2 model or two DM (in the Z3 model). Thus, the three types of events SMu, SM12 and 
SM22, i.e., distinguishable from the identity of SM visible particles, have different edges since they 
correspond to the E'2, E'^, and £'4-type events in a Z2 model or (as mentioned in previous section) 
E2, E-^, and -B4-type events in a Z3 model. This case can be further subdivided into one and more 
than one visible particles in each decay chain. 

In the case with one visible particle per decay chain, one may distinguish Z2 and Z3 models by 
examining the shape of the above three Mt2 distributions. The idea is that, as explicitly mentioned 
above, E2, E!^, and E'^ all have only two DM/DM-like particles, i.e., two DM, one DM and one 
DM', and two DM', and therefore, they have a similar shape of the Mt2 distribution as the £'2-type 
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event. The implication of this observation is E2, E'^, and E'^ all give a sharp upper edge in the M'i'2 
distribution. On the other hand, in the decay chain with two DM in £"3^ 4-type events of Z3 models, 
in general (i.e., away from edge of Mt2), the two DM are not cohinear so that their invariant mass 
of two DM is not 2m£iM, in fact, this invariant mass is not even fixed. Thus, even if the above 
three Mt2- edges for Z3 models are identical to those for Z2 models with two different DM (with 
the second one being twice as heavy as first one), the shapes are not expected to be similar. In fact, 
£"3 and £4 type events in Z3 models give a (relatively) longer tail as already discussed in Sec. 13.31 
Hence, if one of the Mt2 distributions for SMu^ 12, 22 events - again, corresponding to the three 
different topologies - has a sharp upper edge and two of which have a longer tail, then it is likely 
that such events originate from models. 

On the other hand, once there exists more than one visible particle in each decay chain, the 
shape is no longer a useful discriminator. The reason is that, in general, clear sharp edges in the 
Mt2 distributions are not expected here (unlike the cases with one visible particle per decay chain), 
i.e., the number of events/statistics at the MT2-upper edge is small in this case: in turn, this feature 
is due to more constraints which need be satisfied (see Sec. E]). Instead, we can take the advantage 
of "kink" in the plot of M^2^ versus the trial DM mass, which allows us to determine the masses 
of mother and DM particles separately. Using the SMu events, one can evaluate M and moM as 
mentioned in the previous section, assuming that it is a Z3 model. Then we predict the locations 
of the upper edge and the locations of kink for the SM12, 22 events. If the underlying physics is a 
Z2 model (with two different DM particles) instead, then these predictions do not match with the 
experimental results from the associated Mt2 analysis. This is because, in general, the mass of the 
second DM-like particle, m'j-^j^j is not equal to twice of the DM mass, 2mDM- In other words, the 
cross-checking of mother and DM masses between SMu and SM12, 22 events enables us to separate 
Z3 models from Z2. 

5 Identical Visible Particle(s) in the Two Decay Chains 

Next, we consider the case of the visible particle(s) in the two decay chains with one and two 
DM (for Z^ models) being the same. In this case, in the Z3 models, we can not separate £2,3,4 
type events using simply the identities of the visible particles, i.e., SMi = SM2 in Fig. [2] (unlike 
in the previous section). Obviously, we add the three (i.e., E2, £3, and £4-type) distributions of 
Mt2, whose behaviors were discussed above (for non-identical case), to obtain the observable Mt2 
distribution in Z3 models. Of course, for Z2 models (which could potentially fake the Z3 models), 
there are then only £2-type and possibly £2 + events that we discussed earlier. 

If we have only one visible particle in each decay chain, the Mt2 distribution for the £3 type 
events always (i.e., for all m, cf. more than one visible case discussed below) has a lower M™^^^ than 
for the £2 type events (see the left panel of Fig. [3D, so that in principle their addition/combination 
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would give rise to a "kink" in the Mt2 distributioiQ (again for model, but not for Z2 model) 
It turns out, however, that the visibility of this kink is not clear because the Mt2 distribution 
for £"3 type events has a longer tail (that for i?2-type events) as discussed in Sec. [3] (see the right 
panel of Fig. [5]) . It also turns out that the kink will get further smeared out once uncertainties in 
measurements are taken into account. In other words, this kink is not evidence for Z3 models since 
it could be faked by statistical fluctuations in the distribution or experimental errors. 

On the other hand, if there exists more than one visible particle in each decay chain, the Mt2 
distribution for the Z3 model shows somewhat different behavior. As discussed in detail in Sec. 14. 2| 
for trial DM mass above the true DM mass, the upper edge of Mt2 distribution for E2 and E^, type 
events is the same. However, the upper edge of Mt2 distribution for type events is increasingly 
lower than that for E2 type for trial DM mass below the true DM masq^. We therefore expect a 
"moving" kink (as we vary m) in the Mt2 distribution ~ such a kink starts to appear for trial mass 
below the true DM mass (i.e., no kink would appear in the total Mt2 distribution for larger trial 
masses) and the gap between the kink position (i.e., corresponding to the E^ edge) and the overall 
upper (i.e., E2) edge is increasing as the trial mass becomes smalleio. This "moving" feature of the 
kink in the total Mt2 distribution can be further (i.e., beyond simply existence of kink) evidence 
for the existence of E^ type events and thus a proof of Z3 models. However, even though the kink 
is "moving" , it is still hard to identify it in the M'T2 distribution, and thus we do not rely on these 
kinks as a way to distinguish between Z2 and Z3 models. 

These observations motivate us to introduce new methods to separate E^ type events from E2 
type events0 In the following (two) subsections, we develop such a method, and then we apply 
them to the two specific cases (i.e., one visible particle and more than one visible particles in each 
decay chain), and see how to use them to distinguish between Z2 and Z^ models. 

5.1 Separating E2 and £"3 Type Events using Pt/Ht Ratio 

To separate E2 and E^ types of events, we can utilize the fact that E2 type events have one DM 
per decay chain, and £'3 type events have one DM in one decay chain and two DMs in the other 
decay chain. In other words, for E^ type events, the visible particle(s) in the decay chain having 
two DMs in the final state carry less momentum/energy (than in the other decay chain). Thus 
the ratio between the momentum/energy of visible particle(s) on the two decay chains is expected 
to be (relatively) sizeable on average (compared to E2 type events). In order to find out how this 
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■"Such a kink in the Mt2 distribution is not to be confused with that in the plot of M™^"' as a function of m. 
Adding events of £4 type, i.e., two DMs in each decay chain, will introduce another, but even less visible, kink. 
^■^To be more precise, the gap between the two edges relative to mother/DM masses increases. 
^■^In the one visible particle case, the gap between the edges in E2. 3-type events, again relative to mother/DM 
masses, is roughly constant with trial mass so that kink in Mt2 distribution does not move. 

^*Once we separate these two types, we can repeat the program described in detail in Sec. [l] i.e., either determine 
the masses of mother and DM particles from the the upper edges of A{t2 in these two types of events for the case of 
one visible particle in each decay chain or predict the upper edge for Es-type events (using measurements in i52-type 
events) for the more than one visible particle case. 



22 





2 


3 


4 


5 


6 


7 


8 


9 


10 


E2 


0.3375 


0.1629 


0.0957 


0.0612 


0.0413 


0.0285 


0.0225 


0.0178 


0.0144 


E3 


0.7929 


0.6105 


0.4649 


0.3696 


0.303 


0.2525 


0.21 


0.1774 


0.1522 


E2 + neutrino 


0.5366 


0.3323 


0.2319 


0.1733 


0.1342 


0.1101 


0.0887 


0.0738 


0.0628 



Table 1: The percentage of surviving events in E2, and E2 + 1' events for different choice of -Rp™ 
for the case with one visible particle per decay chain. The mother mass is M = 400 GeV and the 
DM mass is rrioM - 150 GeV. 

intuition plays out in real situations, we begin with the case where there exists only a single visible 
particle per decay chain, and then move on to the case where there exists more than one visible 
particle per decay chain. 

5.1.1 One Visible/SM Particle in each Decay Chain 

For the case with one visible particle per decay chain, we consider the Pt ratio of the two visible 
particles as follows: 

^max 

Rpt = (33) 

where p™^^ is the larger Pt of the two visibles coming from two separate decay chains, and p™'" 
is the smaller one. From our physical intuition mentioned above, we expect the Rp^ for E3 type 
events to be larger (on average) than that for E2 type event J^. To verify this expectation, we did 
a simulation using MadGraph/MadEvent for a (toy) model with mother mass M = 400 GeV and 
DM mass moM - 150 GeV. The Rp^ distributions for E2 and E^ type events are shown in Fig. [6j 
We can see clearly that Rp^ for £"3 type events is generally larger than that for E2 type events. For 
comparison, we also included the Rp^ distribution for E2 + V events (i.e., two DM particles and an 
extra invisible, but massless, particle in the full event) in Fig. O One can easily see that the E2 + 1' 
events have on average larger Rp^ than that for pure E2 type events, but smaller Rp^ than that for 
i?3 type events. This observation agrees with our expectation: in the E2 + f events, we have only 
one extra massless invisible particle (relative to £'2-type events), while in E^ type events, we have 
one extra massive invisible particle so that the disparity between the visible particle momenta on 
the two sides in the former case should be relatively smaller. 

Because of different Rp^ distributions of E2 and E^, type events, we can try to distinguish them 
by doing a cut Rp^ > i2p™. The percentage of "surviving" events in E2 and E^ type events according 
to different choice of i?™" is shown in Table [TJ For comparison, we also include E2 + u events in 
Table [TJ We can see that the survival rates for E2 type events are fairly independent of mother 
and DM masses, since in this type of event the energies of the visible particles in two decay chains 
are always comparable. Therefore, 



^We neglect _E4-type events and comment on this issue later. 
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Figure 6: Rp^ distributions for E2 type events (upper-left panel), £^3 type events (upper-right panel) 
and E2 + V events (lower panel) for the case with one visible particle on each decay chain. The 
mother mass is M = 400 GeV; the DM mass is mpiM - 150 GeV. 



• if the survival rates for the events (after the Rp^ cut) are much larger than that of E2 values 
shown in Table [H then we can conclude that the events are not purely E2 type, i.e., it is an 
evidence for existence of another /third invisible particle (whether massless or massive). 

In general, the survival rates for E^ type events are much larger than those of E2 type events, 
but the survival rates of E-^ type events depend on the mother and DM masses. In addition, the 
survival rates for E2 + 1' events are larger than that of E2 type events as well (even though they are 
generally smaller than that of £^3 type events). In this sense, an observation of large survival rates 
might not (by itself) provide a strong support that there exist E^ type events in the sample. 

To get further confirmation of E^-type events (and thus to distinguish Z2 and Z3 models), we 
can employ the Rp^ cut as above and then study the Mt2 distribution of the surviving events. The 
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key idea is to compare the upper edges of the Mt2 distributions before and after the Rp^ cut. If 
the underlying physics model is Z2 type, then clearly we can only obtain E2 type events (or E2 + f) 
events before and after the cut, and the upper edge of its Mt2 distribution is not altered. However, 
if the underlying physics model is Z3, then (before the cut) the total events are an admixture of 
E2 and £'3 type. Since the upper edge of (and i?4)-type events is smaller than those of £'2-type 
events, the upper edge of Mt2 distribution (again before the Rp^ cut) should be that of E2 type 
events. On the other hand, after the cut the upper edge of distribution should be lower than 
before since the surviving events are mostly of £"3 type. 

To illustrate this technique, we apply the analysis outlined above to the previously simulated 
events using a model with M = 400 GeV, mpiM - 150 GeV, and we pick the trial mass to be m = 25 
GeV for the purpose of illustration. Based on the survival rates shown in Table [H we choose 
-Rp™ = 5. Before we do the analysis, we need to investigate whether the Rp^ > 5 cut is "biased," 
i.e., does the cut tend to remove more events with a high Mt2 value For this we consult Fig. [71 
which shows separately the Mt2 distributions for both pure E2 and £3 type events before and after 
the i?Pj cut. By comparing the left panels and right panels in Fig. [7] we can easily see that the 
upper edges for both E2 and E3 type events do not get modified after the Rp^ cut, which suggests 
that the Rp^ cut is not "biased", o In addition, we see that these upper edges in simulated events 
approximately agree with the theoretical predictions (shown by vertical lines)@. 

Now we can demonstrate how to distinguish Z2 and Z3 models using a combination of Rp^ cut 
and Mt2 distributions. First we consider a Z3 model, where we assume that the branching ratios 
for the mother to decay into one DM (and visible particle) and into two DMs (and visible particle) 
are both 50%@ Since we assumed that the visible particles in these two decays are identical, we 
have to combine the M'j'2 distributions for E2 and £3 type events in a 1 : 2 ratio to get the total 
Mt2 distribution. The result is shown in the left panel of Fig. [HI As expected, we can see from 
this figure that the combined E2 and £3 type events have an upper edge in Mt2 distribution that 
agrees with the theoretical expectation for E2 type events (the red solid line). As discussed earlier, 
we can also see that the total Mt2 distribution has a kink near the theoretical M^2^ for £3 type 
events, but it is hard to identify such a kink because of statistical fluctuations. The right panel 
of Fig. [Sj shows the Mt2 distribution for the combined events after the Rp^ cut. It can be seen 
clearly that the upper edge of Mt2 distribution gets reduced. In fact, the new edge agrees with 
the theoretical expectation of the MJ^2^ of £3 type events (the blue solid line). This confirms our 
expectation that the Rp^ cut discards most of E2 type events while retaining a sizeable fraction of 



■^^If the answer is affirmative, then M™J"' after the cut will be reduced even for purely E2-type events. 
^^But a choice of larger 7?p'" would introduce bias in the cut. 
But note that in experiments these predictions are a priori unknown since we do not know the masses of the 
mother and DM separately. 

^^In general, the decay into two DM should be phase-space suppressed relative to the decay into one DM. However, 
in some specific models, this suppression (for the decay chain with two DM) could be compensated by larger effective 
couplings in that chain so that the two decay processes can have comparable branching ratio as assumed here. 
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E2 events before the «p, cut '•^t"" «p> ™' 




AfT2 ^T2 

Figure 7: distribution for £'2 (top two panels) and £"3 type events (lower two panels) for 

simulated events using a model with M = 400 GeV and w-dm = 150 GeV and one visible particle per 
decay chain. The trial mass is chosen to be m = 25 GeV. The left panels are the Mt2 distributions 
before the Rp^ cut, and the right panels are the Mt2 distributions after the Rp^ cut. The solid red 
(dashed blue) lines are the theoretical prediction for the upper edges of Mt2 distribution for E2 
and E's type events, respectively. 



E3 type events, i.e., the events which pass the Rp^ cut are mostly E3 type. 

Of course, we do not know a priori where the M^2^ for £'3 type events (solid blue line) lies due to 
the lack of knowledge of the mother and DM masses. Rather the idea is that we can simply measure 
the upper edge in Mt2 distributions after the Rp^ cut (again, this approximately corresponds to 
that of £'3-type events). Combining this edge with that before the cut (i.e., corresponding to 
i?2-type event) then allows us to evaluate the masses of mother and DM particles as described in 
detail in section I4.]|p| . A complication arises (as follows) in obtaining the M^2^ for events after 
the i?Pj cut. As can be seen in the right panel of Fig. [HI there are still some events beyond the 
theoretical value of M^2^ for £3 type events (the blue solid line in the plot), for example from (a 

^°Note that this cannot be done in Z2 models as discussed in section 14.11 
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Figure 8: Mt2 distribution for combined E2 and type events (1:2 ratio) before (left panel) 
and after (right panel) the Rp^ > 5 cut for the case with one visible particle per decay chain. The 
mother mass is M = 400 GeV and the DM mass moM - 150 GeV. The trial mass is chosen to be 
m = 25 GeV. The solid red (dashed blue) lines represent the theoretical predictions for the upper 
edges of Mt2 distributions for E2 and -Es-type events, respectively. 

small number of) £'2-type events which passed the cut. So we need an algorithm to get rid of those 
"background" events and do a fit to the Mt2 distribution in order to find AI^2^. The details of 
the method we employed are discussed in App. [Cj We apply the above techniques to the simulated 
events. The values of mother and DM masses we obtained from this analysis are 394 ± 8 GeV and 
142 ± 13 GeV, which agree quite well with the theoretical values (400 GeV and 150 GeV). However, 
we expect that uncertainty in energy measurements would introduce additional errors so that a 
more thorough analysis taking into account these effects (which is beyond the scope of this paper) 
is needed in order to be more realistic. 

For comparison, we consider now Z2 models. In these models, we have either E2 ox E2 + V events. 
For pure E2 type events, we have already shown in the upper panels of Fig. [7] that the upper edge 
of Mt2 distribution is not reduced after the Rp^ cut (note that this would not be true if the Rp^ cut 
is "biased"). For completeness, we also consider Z2 models where the mother can decay into one 
DM or one DM plus neutrino, with the visible particle in the two decay chains being identical. We 
again assume that both branching ratios are 50%. Thus, we will obtain a combination of E2 and 
E2 + z^-type events. The Mt2 distributions in this case are shown in Fig. [9l The left panel shows 
that before the Rp^ cut, the upper edge of Mt2 distribution agrees with the theoretical prediction 
of E2 type events and does not have a kink. And, by comparing with the right panels of this figure 
(i.e., after the Rp^ cut, when mostly E2 + i^-type events survive) we can see that the location of the 
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Combined Ei and E2+V events before the ifp, cut Combined £2 and E2+v events after the Rp, cut 
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Figure 9: Mt2 distribution for combined E2 and E2 + v events (1 : 2 ratio) before (left panel) 
and after (right panel) the Rp^ > 5 cut for the case with one visible particle per decay chain. The 
mother mass is M = 400 GeV and the DM mass mpiM - 150 GeV. The trial mass is chosen to be 
m = 25 GeV. The solid red (dashed blue) lines represent the theoretical predictions for the upper 
edges of Mt2 distributions for E2 and E's-type events, respectively. 



upper edge for Mt2 distribution also does not change. These two observations are easily explained 
by the fact that, as discussed earlier, the M^2^ for E2 + z/-type events is the same as for (purely) 
£'2-type events. Based on the above discussions for Z2 and Z3 models, we conclude that 

• by observing whether the upper edge of Mt2 distribution changes (in particular, reduces) 
after the Rp^ cut, we can distinguish between Z2 (including those with neutrino appearing in 
decay of a mother) and Z3 models. 

5.1.2 More than One Visible/SM Particles in Each decay Chain 

Next let us consider the case with more than one visible/SM particles per decay chaiij^ To be 
specific, we consider the case with two visible particles per decay chain. A similar analysis can be 
done if there are more than two visible particles. To separate E2 and E^ type events, we consider 
the ratio of Ht, where HI - ^„ {Pt^l is the scalar sum of Pt^s of visible particles in the same decay 
chain (assuming we know which particles come from which decay chain), and i - I, 2 is the index 
for the decay chains. Ht gives a measure of how energetic the visible particles are in each decay 

•^^Note that in this case, we will get a double-edge in the visible invariant mass distribution from a single mother 
decay, which can already be used to distinguish Z3 from Z2 models [36) . However, it is always useful to have more 
techniques - such as the one, using decays of both mothers in the event, that we are developing here - for such 
discrimination. 
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E2 


0.1542 


0.0401 


0.0142 


0.0065 


0.0028 


0.0015 


0.0009 


0.0006 


0.0005 


E3 


0.7990 


0.5181 


0.3103 


0.1905 


0.1228 


0.0838 


0.0583 


0.0422 


0.0314 


E2 + neutrino 


0.3086 


0.1192 


0.0551 


0.0279 


0.0166 


0.0104 


0.0074 


0.0052 


0.0035 



Table 2: The percentage of survival events in E2, E3 type events and E2 + v events for different 
choice of -R^™ for the case with two visible particles per decay chain. The mother mass is M = 400 
GeV and the DM mass is mo^i - 150 GeV. 

chain. We define the Ht ratio as follows 

trmax 

Rh, = (34) 

where H™^^ - max{Hl,Hf) and i^™"- = m.in{H^ , Hf). From similar reasons to the one visible 
particle case discussed above, we expect Rnt for E3 type events to be larger than that for E2 type 
events on average. To illustrate this feature, we simulate E2, E^ and E2 + v events for a model 
with M = 400 GeV and moM = 150 GeV using MadGraph/MadEvent. The results for the 
distribution for different types of events are shown in Fig. [TOl It can be seen that these distributions 
are very similar to the Rp^ distributions for the one visible particle case shown in Fig. [6l The Rut 
for £'3 type events is on average larger than that ol E2 + y events, which in turn is on average larger 
than that for E2 type events. The survival rates for £'2, £3 and E2 + y events for the cut Rh^ > R^hT 
with different choices of are shown in Table [2j As in the one visible case discussed before, if 
the survival rates for the observed events are much larger than the E2 value shown in Tabled then 
we can conclude that the events cannot be purely E2 type. But large survival rates alone cannot 
be used as a convincing evidence for Z3 models. 

Just like in the case with one visible particle per decay chain, we again use a combined analysis of 
Rhi cut and upper edges of Mt2 distributions to distinguish between Z2 and Z3 models. However, 
there is one major difference between the case with two visible particles per decay chain and the 
case with one visible particle per decay chain. In the latter case, we cannot find the mother and 
DM mass just based on the Mt2 upper edges for events before Rp^ cut. On the other hand, in the 
case at hand, there is a kink structure in the M^2^ vs m plot for £'2-type events which tells us 
both the mother and the DM masses (see section [2]) . And M^2^ for £3 type events is smaller than 
that of E2 type events for fh < m-QM and the same for rh > (see section [321) ■ Thus, before the 

RHt cut, one expects that M^2^ for the combined events is always (i.e., irrespective of the trial DM 
mass) given by that of E2-iype events (just like for the case of one visible particle per decay chain 
discussed earlier). Therefore, in the present scenario, we can find out the mother and DM masses 
before we do any RHt cut. We can then predict the edge in AIt2 for £'3-type events, i.e., after the 
RHt cut (again, the events surviving the cut will be mostly £'3-type). 

We now demonstrate an application of the general strategy outlined above. Based on the survival 
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Figure 10: Rnt distributions for E2 type events (upper-left panel), £"3 type events (upper-right 
panel) and E2 + y events (lower panel) for the case with two visible particles on each decay chain. 
The mother mass is 400 GeV, the DM mass is 150 GeV. 
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rates shown in Table [21 we choose R^^^ = 3 in this case. Fig. [TT] shows the Mt2 distributions for 
the simulated pure E2 and type events before and after the Rnt cut. By comparing the left 
and the right panels in Fig. [TTl we can see that the cut does not alter the upper edge of 
Mt2 distribution for both E2 and E-^ type events. Therefore the Rut cut is not "biased" [f^. We 
then consider a Z3 model where the branching ratios of mother decaying into two DMs (plus two 
visible particles) and into one DM (plus two visible particles) are both 50%. We show the Mt2 
distributions for the combined events before (left panel) and after (right panel) the Rnt > 3 cut in 
Fig. [121 As expected, we see that before the Rnt cut, the upper edge of Mt2 distribution agrees 
with the theoretical prediction of E2 type events (shown by the red line). And, the upper edge 
for the Mt2 distribution gets reduced after the Rnt culjffl, which can serve as an evidence for Z3 
model (cf. discussion of Z2 model below). In addition, as mentioned earlier, knowing the mother 
and DM masses from the kink in the plot of Mr^2^ before the cut as a function of trial DM mass 
, we can predict (shown by the blue line) the upper edge for the Mt2 distribution for the events 
that passed the Rnt cut (it is just the theoretical M^2^ for E^ type events). From the right panel 
of Fig. [121 we see that 

• the observed M^2^ for events that passed the Rnt cut does agree with the prediction (cf. one 
visible particle case above), thus providing additional evidence that the underlying physics 
model is Z3. 

We can compare the above result with Z2 models. If we just have pure E2 type events, then the 
RHt cut does not change the upper edge of Mt2 distribution, as already seen in the upper panels 
of Fig. [TTl As in the one visible particle case, we can also consider the case where there are E2 + V 
events in addition to E2 events: we assume that the branching ratios for mother to decay into one 
DM plus neutrino (plus two visible particles) and into one DM (plus two visible particles) are both 
50%. We show the Mt2 distribution before (left panel) and after (right panel) the Rut cut for 
this case in Fig. [131 from which we can see that the upper edge of Mt2 distribution before the cut 
agrees with the theoretical prediction for £'2-type events and that it again does not change after 
the cut (as expected: see similar discussion for the one visible particle case). 

5.2 A Summary of the Analysis and its Limitations 

Now we summarize the analysis needed to be carried out to distinguish between Z2 and Z3 models 
when the visible particles on each decay chain are identical. 
For the case with one visible particle per decay chain: 

^■^However, a choice of higher give rise to bias. 

■^''Here, we have chosen 9 GeV as the trial DM mass, i.e., (much) smaller than the actual DM mass so that 
the difference between upper edges of Mt2 for E2 and £3 type events, i.e., the reduction of the edge after cut, is 
"exaggerated": see section l3^ 
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£, events before the «„, cut ^2 events after the Rm cut 




Figure 11: Mt2 distribution for E2 (top two panels) and E3 type events (lower two panels) for 
simulated events using a model with M = 400 GeV and moM - 150 GeV and two visible particles 
per decay chain. The trial mass is chosen to be ttt, = 9 GeV. The left panels are the Mt2 distributions 
before the > 3 cut, and the right panels are the Mt2 distributions after the Rnt > 3 cut. The 
solid red (dashed blue) lines are the theoretical prediction for the upper edges of Mt2 distribution 
for E2 and £'3 type events. 
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Figure 12: Mt2 distribution for combined E2 and £^3 type events (1 : 2 ratio) before (left panel) 
and after (right panel) the Rnt > 3 cut for the case with two visible particles per decay chain. The 
mother mass is M = 400 GeV and the DM mass moM - 150 GeV. The trial mass is chosen to be 
m = 9 GeV. The solid red (dashed blue) lines represent the theoretical predictions for the upper 
edges of Mt2 distributions for E2 and -Es-type events, respectively. 



Combined £2 and £2+v events before the Khi cut Combined £2 and £2+v events after the Khi cut 




Figure 13: Mt2 distribution for combined E2 and E2 + v events (1 : 2 ratio) before (left panel) 
and after (right panel) the > 3 cut for the case with two visible particles per decay chain. The 
mother mass is M = 400 GeV and the DM mass moM - 150 GeV. The trial mass is chosen to be 
m = 9 GeV. The solid red (dashed blue) lines represent the theoretical predictions for the upper 
edges of Mt2 distributions for E2 and -Ea-type events, respectively. 
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• We first find M^2^ with different trial DM masses (m) for all the events. We can then 
substitute this value into ^JCe2 + \/ + = M^2^ to find the parameter C E2 (see the 
details in section HTTj) . 

• We apply the cut Rp^ > 5 and find M^2^ with different trial masses (m) for the events that 
passed the cut. If we observe that the M^^^ is reduced after the Rp^ cut (as compared to 
before), then we can conclude that the underlying physics model is Z3 type, otherwise it is 
Z2 type. On the other hand, if the M^2^ is not changed after the Rp^ cut, then we conclude 
that the underlying physics model is Z2 type. 

• If we confirmed (as above) that the physics model is Z3, we can then substitute the M^2^ 
(for various trial DM masses) for events after the Rpf cut into \JCez + \/ Ce^ + fh? - M^^^ 
to find the parameter Ce^- Based on the values of Ce2 ^'^^ Ce^, we can find the mother and 
DM masses simultaneously (again see the details in section I^Jj) . 

For the case with two visible particles per decay chain: 

• We first find M^2^ with different trial masses (m) for all events. We then draw a M^2^ 
versus m plot and find the location of the kink. This can give us both the mother and DM 
masses. 

• We calculate the theoretical predictions of M^2^ for type events using the mother and 
DM masses found in the first step. 

• We apply the cut Rjjt > 3 and find M^f^ with different trial masses (m) for events that 
passed the Rnt cut. If the edge in Mt2 reduces as a result of the cut, then we conclude that 
the underlying physics model is Z3. Otherwise, it is a Z2 model. Furthermore, if the new 
M^2^ agrees with the theoretical prediction for £"3 type events found in the second step, then 
we have additional evidence that it is a model. 

In the above analysis, we have ignored the E4 type events. However, including these events 
would not affect our analysis. Specifically, the M^2^ for E4 type events are always smaller than 
that of E2 and E3 type events (see Eqs. ([T2|) and [20]) so that they would not affect the upper edges 
of Mt2 distribution for events both before and after the RpjRnt cut. However, the survival rates 
for events after the cuts might be modified. In any case, we did not use the survival rates alone to 
distinguish between Z2 and Z3 models. 

The above method of separating E2 and E-^ type events using the Rp^ or Rh^ cut has its 
limitations. If the DM mass is very light compared to the mother mass, then the emitted extra DM 
might not carry away as much energy. Thus, in £'3-type events, the visible particles in the decay 
chain with two DM particles can be closer (relative to the heavy DM case) in energy to those in 
the other decay chain. In fact, the DM becomes similar to a neutrino in this case so that the Rp^ 
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or Rh^ distributions for E3 type events should be similar to those for E2 + z^-type events, in turn, 
not much different from that of E2 type events (cf. heavy DM case), and the distinguishing power 
of the Rp^ or Rjjt cut is reduced. Therefore, in order for the Rp^ or Rnt cut to efficiently separate 
E2 and E^ type events (and hence to distinguish between and Z2 models), we need the mass 
ratio between DM {itiom) mother (M) masses to be sizeablj^. 

Finally, we note that the cut on the ratio of momentum/energy on the two sides of the full 
event can also be used - either by itself or in conjunction with edges in Mt2 ~ for the non-identical 
visible particles case (discussed in the previous section) in order to distinguish E2 and -Bs-type 
events. Of course, in that case, just the identity of the visible particles was enough to separate E^ 
from £'2-type events. 

5.3 Signal Fakes by an (Effective) 2nd DM Particle 

Next, we discuss the strategy to distinguish Z2 models with two different DM particles) from Z3 
models, similar to the discussion in Sec. 14.31 However, now we consider the two decay chains with 
one and two DM, respectively, in Z3 models or with the two different DM particles in Z2 models 
having identical visible particles (unlike in Sec. 14. 3p . In this case, there is a modification from 
the case discussed earlier: we obtain one Mt2 distribution by simply adding Mt2 distributions for 
E'2 3 4-type events in Z2 models (and similarly, £'2, 3, 4-type events in Z3 models) . Let us consider 
first the case with a single visible particle per decay chain. As mentioned before, E2, E'^, and 
E'^-type events, i.e., "sub" -distributions in the case of Z2 models, all give a sharp upper edge in the 
Mt2 distribution. This observation leads to the expectation of two sharp "kinks" ~ at the location 
of the "would-be" (smaller) edges of E!^, and -E4-type events - in the middle of the combined Mt2 
distributior@, in addition to the overall upper edge resulting from E2 type event sH. 

Note that we had a similar discussion for Z3 models in the beginning of section [5j However, in 
the case of Z3 models only E2 type events give a sharp upper edge and the other two type events, 
i.e., E3 and E4 give relatively longer tails (albeit with smaller endpoints than £^2-type events) 
so that two kinks in the combined Mt2 distribution from £'3 and E4 type events are not clear. 
Therefore, clear sharp kink(s) in the Mt2 distribution would suggest that the events are the result 
of a Z2 model (as discussed earlier). 

For the case with more than one visible particle in each decay chain, the above idea of using 
kinks in Mt2 distribution might fail since the edge of the sub-distributions is not sharp, even in 
the case of Z2 models (as discussed in section [^?3]) . Instead, we can do cross-checks like in the case 
of non-identical visible particle(s) in the two decay chains (discussed in section IT3|) i.e., we first 
measure the masses of mother and DM particles by examining the location of the kink present in 



^''Of course, we need vtlomIM < 0.5 in order for the decay chain with two DM to be kinematically allowed, 
^^not to be confused with kink in edge of Mt2 as a function of trial mass! 

^^The situation is similar to the double edge signal in the case of off-shell intermediate particles studied in the 



reference [SS] 
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the maximum Mt2 as a function of the trial DM mass for E2 or E'g'type events, and then predict 
the location of edges in the other types of events. Of course, in order to follow this strategy in 
the present case, one must first separate the events which are mixed, i.e., combination of -©2,3,4 
in the case of Z3 models and -©234 ^2 models, into each individual type by applying a Pt/Ht 
ratio cut (as explained in detail in Sec. 15.11 for Z3 models). Note that even mixed events in Z2 
models with a second DM-like particle can be separated by a Pt/Ht ratio cut because the -Bg-type 
events also have an imbalance in the energy /momentum of the visible on the two decay chains 
due to the difference between nijjM and m^^jQ Of course, we can also do a similar separation 
for the case of one visible particle in each decay chain (which was just discussed above) and then 
repeat the strategy which we discussed in section for the case of visible particle in the two decay 
chains being non- identical, i.e., consider the shape of the separated Mt2 distributions in order to 
distinguish Z3 from Z2 models. 

6 Future Considerations 
6.1 On-Shell Decay Processes 

In the previous sections we have focused only on decay chains where the intermediate particles are 
off-shell. We have used Ht cuts to reduce the inherent Z2-like background (events with two DM 
candidates in the final state) to E^-type events (events with three DM candidates in the final state) 
for the case when each of the pair-produced mothers decay to the same SM final state. As a re- 
minder £"3 events do not occur in models with a Z2 stabilization symmetry. An important question 
is whether this method can be applied to decays where the intermediate particles are on-shell. 

For a quick example on how the decay of the mother particles into intermediate particles that are 
on-shell can complicate the analysis of the previous section, consider the pair production of mother 
particles which are charged under a Z2 stabilization symmetry. We will assume each mother decays 
via on-shell intermediate particles into a dark matter candidate and identical SM particles. Now 
suppose the first mother decays into virtual daughters whose mass differences with each other, the 
mother and dark matter are large. The second mother decays into intermediate daughters whose 
mass differences are all small except for one. It is clear, in this latter case, there will be events 
where the transverse momentum of the visible particles can be small. It is also conceivable, for 
those events, that phase space available to visible particles in the first mother decay is relatively 
large. In these events, on average, the Ht ratio can be large and therefore can fake the Z3 signal. 

To make this more clear, let us go a step further to consider the following explicit example. The 
'^^Obviously, one cannot then use this cut on its own in order to distinguish Z-j, from Z2 models. 
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Figure 14: The production and decay of Q in Eqs. (I35p and (I36p . In the panel on the left the 
intermediate particles generating the decay have a mass difference of 200 GeV. On the panel on the 
right, all of the mass differences are set to 5 GeV with the exception of one which is 590 GeV. (See 
Eqs (|40p through (j42p .) The dark matter mass is the same for both plots and is charged under a 
Z2 symmetry. 



production of two mother particles can be 



pp^QQ, 



(35) 



where Q is a new, heavy particle. Q and Q decay to a final state Q qllx via the following 
on-shell intermediate particles 



Q-> Aq 
Bl 
B^IX 



A' q 
A' ^ B' I 
B'^IX. 



(36) 



Here x is the dark matter particle. A, A', B and B' are intermediate particles. The lowercase 
letters are the other SM particles. The important mass differences are 



6 niQA ^niQ-mA 

5 niAB ^mA-mB 

6 tubx - - 



6 tuqa, ^ruQ- ruA' 
5mA'B' - ITT' A' - fnB' 
SrriB'x - ^B' ~ iTT'x- 



(37) 
(38) 
(39) 



If the mass differences, Stuqa, Stuab, Snisx, are relatively large and the differences, Smq^^,, 
SmA'B'i S^B'x^ small then the Ht ratio may be askew. We show this explicitly in Figure 13 
for 



SruQA - SniAB - Stubx- 200 GeV 
SmQj^i - 5 iJiA'B' - 5 GeV 
SruB'x^ 590 GeV. 



(40) 
(41) 
(42) 
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where rriQ - 700 GeV and - 100 GeV. Ht is plotted above for the small mass and large mass 
separation cases. This argument is general and can be applied to many on-shell decay scenarios. 
For this paper, we considered only off-shell decay scenarios. We leave these questions to future 
work. 

6.2 Invisible Z Decay Background Considerations 

Beyond the on-shell considerations described above, in future work we will also consider accep- 
tance cuts and detector effects consistent with the ATLAS and CMS experiments. An important 
irreducible background to our signal with DM in the final state is the emission of Z bosons that 
decay invisibly. The standard way to mitigate the effect of this background is to look for events 
which radiate off a Z boson decaying to electrons or muons. These Z bosons can be reconstructed 
with standard invariant mass techniques. Knowing the Z boson branching fractions, one can then 
reliably estimate the number of invisible Z boson , we will re- 

quire the signal cross section, after detector and acceptance cuts, to be more than 10% of this SM 
background to preclude statistical fluctuations [35] . 

6.3 Resolving Reconstruction Ambiguities in Multi-Jet Events with Large frp 

Thus far, we have not applied our analysis to many-particle final states with large fJr- Clearly, 
we want to be able to distinguish versus E2 events for those scenarios. The problem is recon- 
struction ambiguities are inherent in, e.g., multi-jet events. Consider the pair production of mother 
particles which decay on-shell into a six-jet final state with two SUSY-like dark matter candidates. 
In order to reconstruct the decay chain, experimentalists must consider 6! = 720 combinations of 
jets! Such ambiguities can prevent accurate reconstruction of events with dark matter candidates, 
not to mention distinguishing scenarios with two versus three dark matter particles in the final 
state. Thus, the question of searching for dark matter in events with many visible final states is 
really the question of finding methods to enable accurate reconstruction. 

Consider the difficult case of decay chains with multi-jet final states and large amounts of miss- 
ing energy. Several groups have proposed methods to increase the efficiency of reconstructing such 
events [371 [S]. These methods require the final state jets to have large pT cuts. Because of the large 
Pt cuts, the regime of phase space where the mother particles are highly boosted is selected. These 
methods are therefore essentially variants of the aptly named "hemisphere method." By definition, 
the hemisphere method considers only events where the mother particles are highly boosted. Thus, 

■^^This technique is common and used often with analysis of estimating the sensitivity ATLAS and CMS to invisible 
higgs decay processes. See for example |41l I42j . Also, as an additional note, it is conceivable that the techniques 
introduced in the previous sections can be applied to separate the signal (DM events) from background processes 
that include invisible Z decays. However, since estimating the invisible Z branching fraction from data as indicated 
above is straightforward and reliable, we focus on that method. 
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final state jets resulting from each mother's decay are correlated in different sides of the detector. 
The correlation allows experimentalists to accurately associate the correct final state jets with the 
right decay chain. There are, however, inherent problems with this method: 

1. New particles produced at the LHC are expected to be heavy and therefore produced domi- 
nantly at threshold. Requiring boosted mother particles pushes production into the regime of 
phase space where the number of signal events are likely to be smaller than that at threshold. 

2. Reconstructing multi-jet events requires hard jets with large cone separations {AR > 0.4 or 
0.5) between the jets. When the jets are correlated in different hemispheres of the detector, 
this requirement will further reduce the number of events which will pass to become signal. 

With the signal events expected to be small (because of the heavy mother production), it is crucial 
to keep as many of the signal events as possible. Thus, the true goal would then be to accu- 
rately reconstruct these events with only the acceptance and kinematical cuts needed to identify 
the signal and reduce the SM background. We now briefly outline [IS] a method to reduce the com- 
binatorics needed to reconstruct events with multi-jet final states with large fJr- We also present 
very preliminary results. The full results will appear in [l5]. To simplify the discussion in the 
next sections, we assume pair production of mother particles with each decaying into identical dark 
matter candidates plus jetjf§. 

6.3.1 Proposed Algorithm 

The reconstruction algorithm is: 

Step 0: Obtain a statistically significant amount of n-jet events (with a large amount of 
fjr) over the SM background. 

Step 1: From these events, first consider only the events that remain in the limit of the 

hemisphere method, i.e., the final state jets have a large pT- The number of 
events retained in this limit must still be statistically significant^'^ over the re- 
maining background. 

Step 2: In the hemisphere method limit, the topology of the signal events can be ob- 
tained. 

To make Step 2 clearer, consider the example of a six-jet final state. A visual representation of 
the topologies is in Figure 15. Again, we assume pair production of mother particles with each 
decaying into identical dark matter candidates plus the SM; also, in the hemisphere limit, the jets 
are required to be well separated. If, in the limit of the hemisphere method, one mother preferen- 
tially decays into five jets and the other mother decays to only one, then the "5-1" topology can 
■^^We relax this assumption in [45] 
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r t c 

Figure 15: All possible six-jet final state topologies with the assumption of pair production of 
mother particles which decay to the SM and identical DM candidates. The 5-l(left), 4-2 (center) 
and 3-3 (right) topologies are shown. 

be identified. As described in the previous section, the identification is due to the mother particles 
being boosted; the resulting decay products are therefore on different hemispheres of the detector 
and can be readily associated with the right decay chain. Another example is possible with one 
mother decaying into two jets and the other decaying into four jets. This produces a "4-2" topology. 
Thus, Step 2 is crucial to identify which topology (or combination of topologiei^) is generated. 

Step 3: In the hemisphere limit, reconstruct the n-jet invariant masses for the 
topology(-ies) determined in Step 2. This is to find kinematic edges. 

It is well known that heavy mother decays into the SM and dark matter generate invariant mass 
"edges." Consider again the case of the 4-2 topology. Reconstructing the four-jet invariant mass 
associated with the mother decaying to four jets plus dark matter generates an edge. This edge 
is roughly at Mmother - -^dark matter • The measured values of these edges will help to increase the 
efficiency of reconstructing when going away from the limit of the hemisphere method. We discuss 
this more detail in the next section. Also, please note, if the decays of a mother are on-shell, even 
more edges are possible than the ones described above. This will further increase the efficiency of 
reconstruction. We show an example of a kinematic edge in the next section. 

So far the above steps have been done in the limit of the hemisphere method. The next step 
is new and is our primary point: 

Step 4: Reconstruct all of the events with the topology(-ies) determined in the limit of 
the hemisphere method. Knowing these topologies serves as a reconstruction 
template to minimize combinatorics for the events which do not have boosted 
mother particles. 

''^It is possible that multiple signal production processes can generate, e.g., 5-1 and 4-2 topologies. If both topologies 
are determined, then reconstructing the full data set requires the assumption of both. More details are included in |45j . 
Also, if e.g., a 4-2 topology is produced, this means events with four and eight jet final states are also produced. By 
focusing on six-jet final states, we do not consider these events. 
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Thus, our point is the use of the hemisphere method to learn the topologies of the signal event. 
This in turn, will reduce the combinatorics needed to reconstruct the full set of signal events 
when only the cuts needed to reduce the SM background are applied. Consider again the example 
of a six-jet final state. If, in the limit of the hemisphere method, one mother is found to decay into 
two jets and the other four jets, the number of ways to associate the jets with the correct decay 
chain is 15. In the case of the 3-3 and 5-1 topologies, it is 20 and 6, respectively. The invariant mass 
edges determined in the hemisphere method limit will further help to reduce the possible ways to 
associate the various jets with each decay chain even further. We describe some specifics briefly in 
the next section. 

6.3.2 Additional Notes on Implementing the Algorithm 

To begin, for each event, the jets are labeled by decreasing pT- Like in the previous sections, the 
scalar sum of the px of several jets, Ht, is defined as 

Ht - t (43) 

a=l 

for a number of jets, j. We use Ht to go to the limit of hemisphere method. To help associate the 
correct jets with the right decay chain, we also construct the transverse sphericity (St) for n-jets, 
where n can be less than the total number of final state jets. The number of jets included in the 
definition depends on the assumed topology for reconstruction. St is defined as 

Sr - ^ (44) 

where Ai^2 sue the eigenvalues of the 2x2 sphericity tensor 

S^j-Zp-^p""' (45) 

K 

where k runs over the number of jets included. In general, when the jets used to compute St are 
collinear or back-to-back, then St 0. Because of this property, St is used to help determine the 
topologies of the signal events in the limit of the hemisphere methocQ. 

In the previous section, we described how reconstructing invariant masses in the hemisphere limit 
can help reduce combinatoric ambiguities. We want to give a basic example of how this is possible. 
Consider again the example of a six-jet final state. Suppose the signal topology is a 5-1 topology. 
In this example, the signal event is generated by pair production of 700 GeV heavy quarks, Qh, 
one of which decays into five jets and 100 GeV scalar dark matter, x. The other mother decays to 
only one jet and the same dark matter candidate, 

pp^ QhQh ^ qq qq qq XX*- (46) 



*'^Beyond, Ht and St other, more common, kinematic variables also will be used identify the topology(-es) in the 
hemisphere limit. We expand on this in detail in |45| . 
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Figure 16: Signal process described in Eq. (jlH]) for a 5-1 topology. The left panel is the five-jet 
invariant mass in the limit of the hemisphere method. Here the hemisphere method is defined 
as the cuts in Eq. ()47p with the exception that all of jets have px > 150 GeV. The right panel 
considers the same events with the cuts in Eq. (j47p . Note: The reconstruction ambiguities obscure 
the kinematic edge. All of the events above the 600 GeV edge are events with poor reconstruction. 

The virtual particles which contribute to the five-jet decay chain have the masses of 580, 460, 340 
and 220 GeV. The five-jet invariant mass has a kinematic edge at approximately Mq^ - - 600 
GeV. In the left panel of Figure 16, we plot the five-jet invariant mass in the limit of the hemisphere 
method. Each dot represents an event. It is clear an edge has appeared at 600 GeV. In the 
right panel, we plot the same five-jet invariant mass with only the cuts needed to reduce the SM 
backgrounds: 

PT leading jet > 100 GcV PT non-leading jets > 50 GeV 

|??jets|<2.5 It > 100 GeV (47) 

5Talljets>0.2 ^>0.2A4fr 
ARjj > 0.4 

where effective mass is defined as 

M^s^Y.PT]cts + fk. (48) 

This cuts are consistent with the ATLAS and CMS collaborations. With only these cuts, we have 
reconstructed the multi-jet events assuming the 5-1 topology determined in the limit of the hemi- 
sphere method. The key point is the edge determined in the hemisphere method limit means all 
of the invariant masses greater than 600 GeV are mis-reconstructions. As described above, the 
amount of mis-reconstructions can be further eliminated if the decay is on-shell and all possible 
kinematic edges are applied. At the partonic level, the reconstruction efficiency for on-shell decays 
for the above example approaches 80% |45j . 
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The next steps to optimize this method is to include QCD backgrounds as well as a variety of 
signal processes with virtual particles different spins. A major fault with Figure 16 is that it was 
generated at the partonic level. The final study [l5] should put the signal processes through pythia 
to get an estimate of hadronization effects. 

7 Conclusions and Outlook 

The LHC experimental program is ramping up. One of its goals is to test the idea of a WIMP as a 
DM candidate, especially one which arises as part of an extension of the SM at the TeV scale. The 
basic idea is that there is a new (global) symmetry under which the SM particles are neutral, but 
some of the new particles are charged. Thus, the lightest of these charged particles is stable and 
under certain circumstances - it certainly has to be electrically and color neutral - can be a DM 
candidate. 

In general, such an extension of the SM also has other (heavier, but still TeV mass) particles 
which are charged under both the new symmetry and under the SM gauge symmetries (espe- 
cially color). So, these DM "partners" (in the sense that they are charged under the symmetry 
which stabilizes DM) can be copiously produced at the LHC. In turn, these DM partners (or 
"mother" particles) will decay into SM particles and DM, resulting in events with missing energy 
and jets/leptons/photons. The goal then would be to piece together the details of the DM model, 
including mass of the DM and the mother particles, by reconstructing the decay chains which lead 
to such events. 

Clearly, the nature of the DM stabilization symmetry plays a crucial role in the above program, 
in part because it determines the number of DM particles in each decay chain of the mother particle: 
for example, in the case where a Z2/parity symmetry stabilizes the DM (typically) only one DM 
appears in decay of a single mother, whereas in a model, both two and one DM are allowed. Of 
course, each event must necessarily involve two such mothers being produced. 

With the above background, in this paper, we showed how to distinguish models with Z3 
DM stabilization symmetry from Z2. Since earlier work [36] only studied measurements of the 
visible part of a single mother decay, here we focussed on using the information from both mother 
decays, including the missing energy in the event which is shared between the two decay chains (or 
sides). We found that the variable Mt2 is useful for this purpose. For simplicity, we studied pair 
production of the same mother, followed by decays to SM particles and DM which involve only 
OjQ^-shell intermediate particles (i.e., which are heavier than their mother particle). Furthermore, we 
assumed that we know the division the SM particles in the full event into two groups corresponding 
to each mother [but of course we do not know (a priori) how many DM are in each decay chain]. 

Clearly, in a model, the events can be classified into three types depending on the total 
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number of DM particles (i.e., two, three or four) vs. only two DM particles for Z2 model. We 
showed that the edges of the Mt 2 distributions are different in these three types of events in a Z3 
model, again even if the same mother is produced (vs. only one edge for Z2 model). This feature 
allowed us to distinguish Z3 from Z2 models. Moreover, we gave predictions for the values of the 
edges in the two new cases, namely, three and four DM in each event, as functions of mother and 
DM mass. Thus, we can extract the mother and DM masses separately using the measurements 
of these different edges for a model. This achievement is especially noteworthy for the case 
of single visible particle in each decay chain since a similar measurement of the mother and DM 
masses is not possible in a Z2 model, based solely on using Mt2 variable!^ 

We emphasized that there are two subcases in the above analysis, namely, the visible particles 
in the decay chain with one DM being identical or different (respectively) to those in the decay 
chain with two DM (for Z^ models). In the case of the visible particles not being identical, it 
is easy to separate the events of the three types mentioned above so that one can then plot the 
respective Mt2 distributions. However, in the case of the visible particles being identical, one 
obviously has only a single Mt2 distribution (i.e., combination of the above three types) to begin 
with. Therefore, we developed a new method to separate out the candidate events with three DM 
vs. two DM in this case, using the observation that the visible particles on the side with two DM 
will have smaller energy/momentum that the visible particles on the side with one DM in the same 
event. This feature is to be compared to the visible energy/momentum being more "balanced" in 
the case of one DM on each side. 

We observed that the above imbalance in the energy/momentum on the two sides hy itself 
provides a hint for the appearance of three DM in the event. However, combining it with edges in 
Mt 2 distributions provides a more powerful discriminator. Finally, we briefly mentioned the case 
of intermediate particles in the decay chain being on-shell. 

Signal Fakes: We pointed out that models with Z2 as a dark matter stabilization symmetry 
can fake the three decay topologies arising in models. An example of such a fake comes from the 
decays in Z2 models which involve effectively a second DM particle (denoted by DM'), for example, 
a heavier, neutral Z2-odd, but collider stable particle. In such a Z2 model, three decay topologies 
are possible if the pair-produced mother can decay into either DM or DM', i.e., we can have DM 
in one side and DM' on the other side, DM on both sides or DM' on both sides. Nevertheless, we 
suggested how to discriminate Z2 and Z^ models by using various strategies. 

Another possibility for a fake of Za-like signals by a Z2 model comes from pair-production of 
more than one type of mother particles which decay into visible particles and, for simplicity, a 
single DM particle. Since the upper edge of Mt2 distributions depends on the mass of the mother 
particle, the introduction of different types of mother particles can easily generate different locations 

''■^The exception to this rule for Z2 model is when the total transverse momentum of the two mother particles is 
non-zero (for example, due to initial/final state radiation). 
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of upper edges so that it could instead be mistaken as the multiple decay topologies of a Z3 model 
(again, with single DM particle). This type of degeneracy between Z2 and Z3 models, however, 
can be easily resolved by taking into account the fact that each event in such a Zi model is still 
"balanced" in the energy/momentum of visible particles on the two sides (again, assuming that 
there is no second DM-like particle) while events in Z3 models can have imbalance if they involve 
one DM on one side and two DM on the other. 

Obviously, this and the earlier work [36] are to be viewed as the first steps in the broad 
program of distinguishing various DM stabilization symmetries using collider data. To this end, 
we outlined some future steps to implement this program. We leave more detailed studies along 
these lines, including relaxing some of the assumptions outlined above and further studying various 
experimental issues, for future work. 
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A The location of M^^2^^ 

In this appendix we will derive the analytic expression for the location of the upper edge in the 
Mt2 distribution. We begin with deriving the general expression of the Mt2 solution for a given 
set of kinematic configuration, then move on to obtaining the maxima of the balanced/unbalanced 
Mt2 solutions, and close with giving the global maximum of the Mt2 distribution, followed by a 
simple application. 

A.l The General Expression for the Mt2 Solution 

The usual Mt2 variable [51 1^ is defined as a generalized transverse mass such that each of pair- 
produced mother particles decays into visible particles and one dark matter particle of the same 
type. However, we do not restrict ourselves to such cases, i.e., we extend our consideration to the 
cases with more than two DM in a full decay chain (e.g., E-^ and type events in Z3 models). 
Nevertheless, in the analysis of Mt2 variable, we still hypothesize that two dark matter particles 
with equal mass (i.e., one DM per chain) are involved in the full decay process, i.e., we employ the 
"naive" Mt2 method (as mentioned at the beginning of section [3|) . 

The left diagram of Fig. [T7] illustrates the decay process of pair-produced mother particles 
that we are taking into consideration. Here the "blob" denotes off-shell intermediate particles or 
an (on-shell) point interaction. The red dashed lines represent any particles charged under dark 
matter stabilization symmetry whereas the black solid lines represent any visible/SM particles. M 
is the mass of the mother particle, which must be charged under DM stabilization symmetry. As 
mentioned above, each mother particle can decay into the multiple number of invisible/DM particles 
as well as the multiple number of visible/SM particles, and this extended possibility is explicitly 
depicted by the multiple number of red dashed and black solid arrows behind the two small blobs. 
Each visible/invisible multi-particle state can be collapsed effectively to a(n) visible/invisible single 
particle state by introducing invariant (transverse) mass, which will be manifest in the detailed 
formulas later. In this sense, and rnf'^ can be understood as the total invariant masses 

formed by visible or invisible particles belonging to the same decay chain. 

For defining the Mt2 variable one should note that we are not aware of the DM mass in advance. 
Hence, the best we can do is to introduce trial DM mass. Since we perform the naive Mt2 analysis 
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Figure 17: The left panel illustrates the decay process of interest which pair-produced mother 
particles go through. M, mi"'\ and m-'*^ (a = 1, 2) denote the mass of mother particle, total 
invariant masses of visible particles and invisible particles in the same decay chain, respectively. 
The right panel illustrates the effective configuration of such a decay process. 



as mentioned above, i.e., we assume a single type of DM in each decay chain even if the actual 
physics could be different, we employ only one type of trial DM mass rh and construct the 
variable as follows [21 [3] : 

MT2fp?'\m;('\p;('\ml('^m)E min Tmax (m|'\ m|'H1 (49) 

Here each transverse mass of the decay product M^"^ (a = 1,2) is given by 

( a4") )' = (m^(") )' + + 2 (i?;^'^) - p^(") . p(?) ) (50) 

where m!^"'\ Py"\ and E^"'^ are the total transverse invariant mass, transverse momentum, and 
transverse energy of visible particles: 

o a>/3 

Pf^ - Ep:? (52) 

a 

E<-^ . (53) 

a 

and and E^,^'*^ are the transverse momentum and energy of the (assumed-to-be-one) trial DM 
particle in each decay chain. Here rn^"'^ indicates the mass of ath visible particle in ath decay 
chain (o = 1, 2) and t^^""^ indicates the (1+2) momentum on the transverse plane which is defined 
as 

,v(a) ( T-\v(a) v(a)\ I ( v(a)\^ I v(a)\'^ v(a)\ /r-A\ 

h =[EiT^'PiT )-\y(viT ) )'P/r j' (54) 
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Trial DM transverse momentum 



Trial DM transverse momentum 



Figure 18: Graphical configurations to give rise to a balanced Mt2 solution (the left panel) and 
an unbalanced Mt2 solution (the right panel). The dashed line indicates the solution for the Mt2 
variable to take. 

and the metric for this type of momentum is diag(l, -1,-1). There arise two noteworthy things: 

• As far as the range is concerned, the transverse and the regular invariant masses have the 
same lower and upper limits. Moreover, since the Mt2 solutions of interest arise at either of 
the two limits, one may consider the Mt2 where m^"'^ are replaced by the regular invariant 
masses of visible particles : 

a a>f3 

• As advertised earlier, the entire visible states in the same decay chain can be understood 
effectively as a single visible particle whose "effective" mass is given by mi"'\ On the other 
hand, the corresponding "effective" quantity for invisible particles m^"^ does not seem to be 
contained in the Mt2 variable. In fact, the Mt2 variable depends implicitly on m^"\ which 
will be cleared shortly. 

From these two observations we can reduce the decay of pair-produced mother particles into two 
multi-particle states (left panel of Fig. [T7j) to an effective kinematic configuration where there exist 
two simple 2-body decay chains shown in the right panel of Fig. [T71 

For such Mt2, there are two types of solution which are called the "balanced" Mt2 solution 
and the "unbalanced" Mt2 solution, and the fact that there always exist some events to give such 
solutions was proven [2l [3l 0] . The balanced solution arises when M^"^ > M^^^ for both a(= 1, 2) 
with the trial DM momentum having the value to accommodate [b t a) at their global (or 
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so-called "unconstrained" ) minimum which will be defined shortly (see the left panel of Fig. [TH]) , 
and otherwise, the unbalanced solution arises (see the right panel of Fig. [T8|)F^ The global minima 
for and are easily evaluated by differentiating Eq. (j50p with respect to the trial DM 

momentum and finding the stationary point [2l [3]: 



«') - 

V /mi 

«') - 

\ /mi 



m^^^ + rh (56) 
m^^^ + m (57) 



The balanced and the unbalanced Mt2 solutions for a given set of the "effective" visible and invisible 
masses which are shown in the right panel of Fig. [T7] are as follows [21 [3lll] : 



/ 

1 + 



(^2_(^(l)„42)j2j ^^^^ 



for 



^J [ 2A-[mi'^y-(mi^^y) 
M^f"' = m + mi^^ (a = 1, 2) (59) 

A^E:}'^EfKp^^^'^-pf\ (60) 



Note that the unbalanced solution is simply given by the unconstrained minimum of M^''^ and for 
a fixed set of visible and invisible masses the balanced solution is bounded above at 

In fact, if we take an adequate number of events, we can always find some event which corresponds 
to such an upper bound [3]. Since we are interested in the M!^ equal to its own upper bound, 
we henceforth assume that A is understood as Eq. (|61|) unless there arises any confusion. Here 
= + (mi"^^ and g"^") is the magnitude of the total momentum of visible particles 

seen in the rest frame of their mother particle. The explicit expression for q""^"-^ can be easily 
determined in terms of the masses of mother, visible, and invisible particles as follows: 

= ^^{(a/ . m<"))'' - (ml-))] {(m - ,n<"')^ - (,n<"')]^ (62) 

One thing to be emphasized is that the dependence of the Mt2 variable on the effective invisible/DM 
mass m^"^^ first appear in qi"'\ In other words, Mt2 is an implicit function of m^'^^ via qi""^ 
as mentioned before. Furthermore, we include the possibility that multiple (massive) invisible 



**If the total invariant masses of visible states in both decay chains are the same, only the balanced solutions arise. 
The reason is because the unconstrained minima for both Mt are identical, there is no possibility that the kinematic 
configuration like the right panel of Fig. [18] is made. As an example, if there exists only one massless visible particle 
in each decay chain, the Mt2 values are always given by the balanced solution. 
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particles are emitted in each decay chain unhke the previous studies (which considered the cases 
with two invisible particles having non-identical masses). Hence rnf'^ as well as have their 
own range once multiple visible and invisible particles are involved in the given decay process. For 
off-shell intermediate particles the respective ranges are given by (See, for example, |39j ) 

^ ml^'^^ . m'^:;>^^ < mi'^) < mi%,^ = M - ^ m^^^^ (63) 

a 13 

The lower limit corresponds to the situation in which particles described by p^^"^ {Pp"^ ) '•'^^ '^^^^ 
in their center of mass frame so that they move with the same velocity in any frame. The upper 
limit corresponds to the situation in which particles described by p'p"'^ {pa"^^ are at rest in the 
overall center of mass frame of the final state described by p'a"'^ and p'p"'^ . 

A. 2 The Maximum Balanced and Unbalanced Mt2 Solutions 

For the decay with visible/invisible multi-particle final states, it is obvious that balanced/unbalanced 
Mt2 solutions have their own range due to the existence of the range in either rri^^ or rnf'^ or 
both of them. As far as the upper edge in the Mt2 distribution is concerned, either the maximum 
balanced or the maximum unbalanced solution appears as the global maximum. For the unbalanced 
solution, one can easily derive the following relationship from Eqs. (|59|) and 

■^Taax.,unbal 



= m + max [m^^^iax: "i^^],^^] . (65) 



For the balanced solution, however, it is not easily seen which values of and rnf^ (a = 
1, 2) will form the maximum balanced solution because of the complication in the corresponding 
expression (^m!^Y given in Eq. (j58p . In order to identify those values, we are required to carefully 
investigate the functional behavior of {Mtj^^^ according to the changes in m-'*^ and mi°'\ which 
will be considered in order. 

A. 2.1 The Change in m^"'^ 

To see the dependence of (mI^^) on m^^\ we simply take the partial derivative: 

d{M^^if D 



(66) 



where 

D = {BC' - B'C)X^ + 2A'VWCX + B'BC 
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with the followmg notations: 

A'^ B'^^^^2A\ C'^^^^2AA'. (69) 

One can easily see that A, B, and C are always positive for any set of mi""^ and m^'^\ and that 
only positive X is allowed by construction. Also, one can easily prove that A' is negative. 
The solutions to -D = are given as follows: 

X, = ^^^^ (70) 

A-(mi'^f 

X2 = ^ 2 - (71) 

(1) (2) / / (2)\'^ 

For TTit, ' ^ m\, ' it can be proven that either A - I m\, ' \ or yl - I m\, ' 1 must be positive and 
the other is positive or negative depending on the parameter space formed by iri^^ and [3]. 
Hence, one of the two solutions given above must be negative, which is unphysical, the other is 
either physically allowed or not. Actually, it turns out that the signs of A-(ml^)) and A- 
are connected to the coefficient of X"^ in D in the following way: 

EC - B'C = 2A' [a - (m,(i))') ( A - (m(2) )') . (72) 

Let us assume that m\, is larger than m„ . In this case, A - [m\, 1 is always positive, i.e., X2 

is always unphysical. Since A' < as mentioned above, if A - is positive as well, then D, 

which is a quadratic function in X, becomes a parabola bounded above, and the two solutions Xi 
and X2 all are negative, i.e., unphysical. Therefore, D < for arbitrary (physically-allowed) X or 
fh. On the other hand, if yl- is negative, then D turns into a parabola bounded below, and 

Xi becomes a physically allowed solution. Therefore, D < for < X < Xi and D > for X > Xi. 
However, in [3j it was shown that Xi gives rise to 

M^""^ {X^Xi)^m + mi^^ , (73) 

which is simply the unbalanced solution for 7n^^^ > ml^\ Moreover, they showed that this implies 
that Xi corresponds to the boundary between the balanced domain and the unbalanced domain. 
In other words, with X being larger than Xi the balanced solution is reduced to the unbalanced 
solution. One can make the same argument and lead to the same conclusion for the opposite 
configuration, i.e., ml^^ < m[P . Also, the dependence on m^'^^ can be easily checked by following 
similar arguments. Based on hitherto observations, we have 

= (..p) = 4^)^,, ) . rnf ) J (75) 
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for any set of mi^^ (a = 1, 2). 



A. 2. 2 The Change in ?nl"^ 



The early work on the dependence of (iVf/^g')^ on m^^'^ was made in [1]. Here we simply provide 
the final results and mention some modification from the original expression. 



< for m < m' 



> for rh>m' 



rh < m' 



("4/^ = mi'^ = m^l^) for 

Mlp^ ijni^^ = m^l^iax, nii^^ = ml^max) for m>m' 



(76) 



(77) 



Here the "kink" location m' can be identified as the true dark matter mass moM if only a single 
type of DM is involved [4]. However, in general, it differs from because we do not restrict 

our consideration to the case with one single-typed DM emitted in each decay chain. Therefore, its 
expression is written in terms of all parameters (i.e., M, m^"\ and m^""^), and it can be calculated 
by solving the following equation [3]: 



2A{A-{m>f)+A-{m<f 

where and denote the heavier and the lighter (invariant) visible masses between the two 
decay sides, respectively, and A is defined as 

A^^^. (79) 
d{m>) 

A. 3 Discussions and Application 

It is a well-known fact that there arises a "kink" in the M™^^^ a function of the trial mass once 
there exists more than one visible particle in each decay chain and its location is at m = m dm for the 
cases with a single identical DM particle per decay chain. For more extended case, i.e., m^^^ ^ m^'^\ 



one can simply solve Eq. (|78|) . In turns out, however, that this is not the only way of obtaining 
the kink location. An alternative and simpler way is to find the intersecting point between the 
maximum balanced and unbalanced solutions. In other words, the solution to satisfy Eq. (j78p also 

then have 

v2 / \2 



satisfies the relation M^^^'^""^ - _ Yot simplicity, let us assume that m^^max > ^^i^^max- We 



=(m + maax) (80) 

and 

^^max.M J2 ^ ^2 ^ ^ ^ ^J ^ (B + Arh^) (81) 
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where A is evaluated at m^""^ = ^-i min discussed before. Letting Eqs. ([80|) and ([8T]) be equated 
and doing some tedious algebra, one can end up with 

^^ = K^fK^)', (82) 

which is valid only with rrii^^ being their maximum. Note that M^^^'^"^^ at rh - m' arises when 
mi""^ - r?T,^,"max from Eq. (j77p . Hence, the above-given relationship holds, and the location of the 
kink can be evaluated by finding the intersection between the maximum balanced and unbalanced 
solutions. 

This observation, actually, leads us to the expressions for M^^^'^"^^ and ]\/[^^^^^'^^°-'- _ Note that 
it was proven that the balanced solution contributes to the upper edge of the Mt2 distribution at 
m < m' in H]. Also, it is straightforward to prove that the maximum unbalanced solution is larger 
than the maximum balanced solution at rh>m' . Therefore, as long as the values to give the upper 
edge of the Mt2 distribution are concerned, it can be (effectively) understood that the maximum 
balanced solutions occur at mi""^ - rni}„- and mS"'^ - mS"'}„- and the maximum unbalanced solutions 

1/ ^ iniii t X , 111! 11 

do at the maximum of the two m^^max - 2) for any fh. As an example, if all visible particles 
are assumed massless, the maximum balanced and unbalanced solutions are given as follows: 



max,ba( 
T2 



M max, unfeai 
j,^ - m + max 



4M2 

?j,max' D.max 



N iM-^ -i^^ -^"^ ^''^ 

[-l'2in, -Sin] (84) 



rh + M - min 



Obviously, the upper edge in the Mt2 distribution is determined by the maximum value among 
many events for a given trial DM mass. 

M^^i'^im) = max [MT2(m)] (85) 
many events 

Based on the above-discussed understanding, one could expect that taking the maximum between 
j^max,bai ^^^^ j^max,unbai ^-j^ result in the Same value as the above-given Eq. ()85p . 



Mmax r 71 j-max,6a/ ii j-max,iin6a/ 

2^2 = max Mrp2 , Mrp2 



(86) 



It turns out, however, it is true only for the case where there exists more than one visible particle 
on each decay chain. In the case where there is only one visible particle per decay chain, one can 
prove that the maximum unbalanced solution is less than the maximum balanced solution for any 
rh so that the M^^^ is simply governed by the M^^'^'^'' ^m^^"^ = —-i mm)- 

B The Existence of a Kink in M^2^ Versus fh 

As discussed in App.fXJ it is obvious that for the cases where there is only a single visible particle in 
each decay chain, the M^2^ as a function of the trial DM mass behaves like a smoothly increasing 
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Figure 19: The kinematic regions to have a kink or no kink in the graph of ^ versus 

curve because the upper edge is solely governed by the "balanced" solution in Eq. ([58]) . However, if 
there exists more than one visible particle per decay chain, the competition between the "balanced" 
and the "unbalanced" solutions, which is explicitly given in Eq. (j86p . gives rise to the possibility of a 
kink (i.e., no longer smooth) in the plot of versus rh. In fact, this approach, the competition 

between the two types of solutions, enables us to examine easily whether or not there exists a 
"kink" on the function of the location of M^^^. 

In order to have a kink in M^2^ as a function of m, the two functions of the maximum balanced 
and unbalanced solutions over the trial DM mass, i.e., M^2^'^'^\rn) and M™2'^^'""''"'(m), must cross 



each other. From Eqs. (j58p and (j59p they are monotonic functions in m, and the slope of Eq. (j58p 
is not greater than that of Eq. (j59p (= 1) over the entire range. These two observations tell us that 
once a crossover is made, no additional crossovers are made. Therefore, it is sufficient to check 
whether or not the relative sizes of their corresponding functional values at m = and m ^ oo 
are flipped for ensuring such a crossover. Let us assume that the visible particles are massless for 



simplicity!^ From Eqs. (j83p and (j84p one can easily prove that larger than ]\{^'^^'^"-^ 



at m -s- oo, and thus M^^^'^"^^ should be larger than ]\^^^^^'^^^°-^ at m = to obtain a kink. Their 
functional values at m = are expressed as follows: 



M2 



M^'^Xm^O) = '-^^ ^ (87) 



M max, wnfoa/ / ~ o\ 
2^2 \fn = 0) = M - mi 

where "m^^,}- = nii and rrif'^-,„ = ?7i2, and we assumed mi < m2 without loss of generality. Therefore, 



'One can easily apply the same argument for the case of massive visible particles. 
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the condition to have a kink is simply given by 



which can be further simphfied to 



M M y M + mi ^ ^ 



Likewise, one can easily find the condition to have no kink as follows: 



^>J^^. (91) 

Fig. [19] shows the regions for a kink or no kink pictorially; the red area represents all possible 
kinematic configurations for the existence of a kink while the gray one all possible kinematic con- 
figurations for no kink. Elsewhere is not physically allowed due to the assumption that mi < 7712. 

There is a special case where each decay chain emits only one single-typed DM particle, i.e., 
mi = m2. The Z2 models or the E2 type events of Z3 models belong to this case. The range to 
satisfy Eq. (|9ip . i.e., the condition to have no kink, is ^ > 1 or ^ < -2 both of which are not 
physically allowed. Therefore, we always obtain a kink in M™J^^ as a function of the trial DM mass 
as expected [1]. 

As another concrete example, let us take E3 type events of Z3 models, where one of the two 
decay chains emits a single dark matter particle whereas the other one emits two dark matter 
particles with intermediate particles ojff -shell. Like before, we assume that all DM particles to be 
emitted in the full decay process have the same mass so that the minimum of the effective dark 
matter mass of the two DM side to give the maximum balanced solution is m^j|^ = Itudm, i-e., 
mi = m£,M and m2 = 2m£,M- From Eqs. (j90p and (j9ip the conditions to have a kink or no kink are 
more simplified to 

A 1 • 1 n 'f^DM - 1 

A kmk: < < (92) 

AI 2 ^ ' 

AT 1 ■ 1 ^3 - 1 ruDM . .„o\ 

No kmk: < < 1, (93) 

2 M ^ ' 

which was mentioned in Sec. 13.21 and demonstrated in Fig. HI 

C Algorithm to Find the Upper Edge of Mt2 Distribution 

In this appendix, we describe an algorithm to identify the M^2^ for events after the Rp^ cut. As 
we discussed in section [5. l.H in Z3 models where there is only one visible particles per decay chain, 
and the visible particles in the decay chains with one DM and two DM are identical, the total 
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Mt2 distribution becomes a combination of the distributions of E2 and £"3 events. So the idea is 
to apply an Rp^ cut (a cut on the ratio of PtS of visible particles on the two decay chains in the 
same event) to "remove" the E2 events. This in principle can give us a relatively pure sample of 
S3 events, which has a smaller M^2^. But in practice/reality, there is still a small number of E2 
events that survive the Rp^ cut. Therefore, the upper edge of Mt2 distribution for events after the 
Rp^ cut is hard to be determined due to the "contamination" of E2 type events. This is shown in 
the right panel of Fig. [8l which shows clearly that there is a small number of events which has 
Mt2 beyond the M^2^ of -E3 type events. Here we propose an algorithm to identify/extract the 
"would-be" M^2^ for i^s events by removing E2 contamination events and then do a fitting to the 
resulting distribution. 

First, we need to "subtract" the contaminating events. To do this, we calculate the moving 
average of the number of events per bin including the last n bins in the Mt2 distribution: An- 
The choice of moving average makes this quantity rather stable under statistical fluctuations as we 
increase n. However, as we increase n to the point below M™^^^ for E'^ type events, we start to get 
a sharp rise on An- Based on this, we define rimax to be the bin such that ^nmax+i ^ 2.5A„^^^. This 
bin is considered as a rough separation point between "contaminating" E2 type events and the start 
of E^ type events. And we treat ^nmax ^ rough estimate for the E2 type events contribution to 
the number of events per bin. 

Next, we pick events with n > nmax and subtract An^^^ from the number of events in each bin. 
This gives us an approximate Mt2 distribution for pure E^ type events. Since we do not have an 
analytical formula for the Mt2 distribution for £'3 type events, we can only do a fitting for events 
near and to the left of the bin nmax to find the upper edge of the Mt2 distribution. We choose 
two fitting functions, one linear function and one quadratic function, and did the fitting separately. 
Our final answer for the M™^^^ is given by the average of the values obtained by these two fitting 
methods, and their difference is regarded as the error due to fitting. 
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